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DNA  methylation  as  an  epigenetic  factor  in  the  development  and  progression  of 
polycythemia  vera 

Jean-Pierre  Issa,  M.D. 

Grant  MP043 15 


Study  Site: 

Department  of  Leukemia,  University  of  Texas  M.  D.  Anderson  Cancer  Center 
1515  Holcombe  Blvd.,  Houston,  TX  77030 

Introduction 

Polycythemia  vera  (PV)  is  the  most  common  myeloproliferative  disorder  with  a  yearly 
incidence  of  28  per  1  million  people  and  a  slightly  higher  prevalence  in  males.1  PV  is 
characterized  by  clonal  expansion  of  erythroid,  myelomonocytic,  and  megakaryocytic 
lineages,  erythrocytosis  being  the  most  prominent  clinical  manifestation  of  PV.2  The 
disease  is  associated  with  a  significant  morbidity  and  mortality,  including  thrombotic 
and/or  hemorrhagic  events,  and  a  risk  of  an  evolution  into  myelofibrosis  and  leukemia.3'5 
An  acquired  activating  V617F  (1849G>T)  mutation  of  JAK2  tyrosine  kinase  has  been 
recently  found  in  the  majority  of  patients  with  polycythemia  vera  (PV),  in  about  half  of 
those  with  essential  thrombocythemia  (ET)  and  myelofibrosis  (MF),6'10  and  in  10-20% 
patients  with  chronic  myelomonocytic  leukemia,  Philadelphia-negative  CML,  atypical  or 
unclassified  myeloproliferative  diseases  (MPD)  and  megakaryocytic  leukemia.11"13  It  is 
not  known  what  other  factors  determine  the  disease  phenotype  of  PV,  MF,  and  other 


W81XWH-05- 1-0535;  PI:  Issa,  Annual  Report  2007,  Page  5 


MPD,  and  what  factors  other  than  JAK2  lead  to  disease  progression.  Very  little  is  known 
about  epigenetic  changes  in  PV.  Epigenetic  lesions  have  been  recognized  to  be  important 
in  cancer,  in  particular  in  older  individuals.  Methylation  of  cytosines  in  the  CpG  sites 
clustered  in  the  gene  promoter  regions  results  in  epigenetic  gene  silencing,  and  acts  as 
one  of  possible  mechanisms  of  tumor  suppressor  inactivation  in  cancer.14  Diverse 
myeloproliferative  phenotypes  caused  by  a  single  point  mutation  of  JAK2  tyrosine 
kinase,  lack  of  other  genetic  specific  lesions  in  PV,  and  its  association  with  higher  age 
lead  us  to  propose  the  hypothesis  that  epigenetic  silencing  may  play  a  role  in  the 
pathogenesis  of  PV. 

Statement  of  Work 

Task  1.  Discover  genes  whose  promoter-associated  CpG  islands  are  methylated  in 
patients  with  polycythemia  vera  (PV),  months  1-18: 

a.  Identify  in  the  M.  D.  Anderson  database  all  patients  with  PV  for  whom 
archived  bone  marrow  biopsies  are  available  (month  1). 

b.  Collect  paraffin-embedded  bone  marrow  biopsies  on  all  patients  (projected 
100  patients,  10  cuts/month,  months  1-10) 

c.  Collect  existing  blood  samples  from  PV  patients  at  M.  D.  Anderson,  and 
from  the  external  collaborator  at  Baylor  College  of  Medicine  (projected 
50-60  patients  per  year,  months  1-36). 

d.  Extract  DNA  from  paraffin  cuts  (start  month  1  -  ongoing  until  all  samples 
collected,  months  1-10)  and  from  blood  samples  (months  1-36). 
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e.  Perform  genome-wide  screening  for  promoter-associated  CpG  islands 
differentially  methylated  in  15  patients  with  polycythemia  vera  in  the 
polycythemic  phase,  15  patients  who  developed  myelofibrosis  and  15 
patients  who  transformed  to  leukemia.  We  will  use  Methylated  CpG 
Island  Amplification  coupled  with  Representative  Difference  Analysis 
(MCA-RDA)  as  a  screening  method  (months  2-18). 

Task  2.  Determine  the  methylation  and  expression  profile  of  candidate  genes  in  the 
polycythemic  phase  of  PV,  patients  who  developed  myelofibrosis  and  patients  who 
transformed  to  leukemia.  Months  2-36. 

a.  Bisulfite  treatment  and  PCR-based  methylation  analysis  for  all  the  genes 
discovered  by  MCA-RDA  and  candidate  genes  involved  in  growth  factor 
signaling  (months  2-20) 

b.  Analyze  samples  for  gene  expression  by  real  time  quantitative  RT-PCR 
(months  13-36) 

c.  Statistical  analysis  of  the  collected  data  (months  21-22) 

d.  Validation  of  the  results  on  prospectively  collected  samples  (months  23- 
36) 

Task  3.  Begin  exploring  the  function  of  the  most  promising  genes  using  in  vitro 
cultures  and/or  transfection  experiments.  Months  13-30. 

a.  Determine  whether  specific  inhibition  of  candidate  PV-methylated  genes 
in  normal  cells  would  mimic  the  PV  phenotype  of  hypersensitivity  of 
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erythroid  progenitors  to  erythropoietin  (months  13-30). 
b.  In  case  the  candidate  genes  are  methylated  and  silenced  in  leukemic  cell 
lines,  we  will  restore  their  expression  using  standard  gene  transfection 
technology.  The  transfected  cell  lines  will  be  examined  for  growth 
characteristics  and  in  vitro  differentiation.  The  effect  of  this  transfection 
on  the  function  of  putative  affected  pathways  will  also  be  examined 
(months  13-30) 


Task  4.  To  assess  the  prognostic  significance  of  aberrant  methylation  in  PV  we  shall 
perform  retrospective  multivariate  analyses  of  the  association  of  CpG  island  methylation 
with  survival  and  probability  of  transfonnation  to  myelofibrosis  or  leukemia  (months  24- 
36). 
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Progress  on  Task  1 


We  proceeded  with  MC A/microarray  analysis  of  PV  samples  to  address  Task  le  as 
planned. 

We  performed  methylated  CpG  island  amplification  (MCA)  in  four  patients  with 
different  types  of  myeloproliferative  disorders:  one  patient  who  developed  PV  from  ET, 
one  patient  who  transformed  from  PV  to  MF,  one  patient  with  ET  transfonned  to  MF  and 
one  patient  with  primary  myelofibrosis.  We  applied  MCA  amplicons  to  high  density 
oligonucleotide  microarrays  (Agilent)  to  detect  genes  that  are  methylated  in  MPD 
patients  and  not  methylated  in  normal  controls.  The  work  and  data  analysis  are  in 
progress.  We  analyzed  the  data  by  Agilent  ChIP  Analytics  1.3.1  software  and  detected 
860  of  10419  (8.25%)  genes  significantly  more  methylated  in  MPD  patients  than  in 
nonnal  controls  (P<0.05).  Top  50  genes  with  the  P  value  of  significance  below  0.01  are 
listed  in  Table  1. 

Table  1,  Genes  differentially  methylated  in  MPD  patients. 


Gene 

Description 

Location 

Type 

B4GALT2 

UDP-Gal:betaGlcNAc  beta  1,4- 
galactosyltransferase, 
polypeptide  2 

Cytoplasm 

enzyme 

CHDH 

choline  dehydrogenase 

Cytoplasm 

enzyme 

ME1 

malic  enzyme  1,  NADP(+)- 
dependent,  cytosolic 

Cytoplasm 

enzyme 

MYH6 

myosin,  heavy  chain  6, 

Cytoplasm 

enzyme 

PDE4C 

phosphodiesterase  4C,  cAMP- 
specific 

Cytoplasm 

enzyme 

SARDH 

sarcosine  dehydrogenase 

Cytoplasm 

enzyme 

HSPB8 

heat  shock  22kDa  protein  8 

Cytoplasm 

kinase 

PRKCZ 

protein  kinase  C,  zeta 

Cytoplasm 

kinase 

ABLIM2 

actin  binding  LIM  protein  family, 
member  2 

Cytoplasm 

other 

LLGL2 

lethal  giant  larvae  homolog  2 

Cytoplasm 

other 

MEGF6 

multiple  EGF-like-domains  6 

Cytoplasm 

other 

PLEKHG5 

pleckstrin  homology  domain 

Cytoplasm 

other 
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containing,  family  G  member  5 

SH3PXD2A 

SH3  and  PX  domains  2A 

Cytoplasm 

other 

TSC2 

tuberous  sclerosis  2 

Cytoplasm 

other 

ACP5 

acid  phosphatase  5,  tartrate 
resistant 

Cytoplasm 

phosphatase 

EGFL8 

EGF-like-domain,  multiple  8 

Extracellular  Space 

other 

EMILIN2 

elastin  microfibril  interfacer  2 

Extracellular  Space 

other 

MBP 

myelin  basic  protein 

Extracellular  Space 

other 

ADARB2 

adenosine  deaminase,  RNA- 
specific,  B2 

Nucleus 

enzyme 

AKAP8 

A  kinase  (PRKA)  anchor  protein 

8 

Nucleus 

other 

NEURL 

neuralized  homolog 
(Drosophila) 

Nucleus 

other 

ZNF541 

zinc  finger  protein  541 

Nucleus 

other 

ZNF646 

zinc  finger  protein  646 

Nucleus 

other 

AIRE 

autoimmune  regulator 

Nucleus 

transcription  regulator 

NFIX 

nuclear  factor  l/X  (CCAAT- 
binding  transcription  factor) 

Nucleus 

transcription  regulator 

TCEA2 

transcription  elongation  factor  A 
(Sll),  2 

Nucleus 

transcription  regulator 

GPR153 

G  protein-coupled  receptor  153 

Plasma  Membrane 

G-protein  coupled  receptor 

TMEM37 

transmembrane  protein  37 

Plasma  Membrane 

ion  channel 

TRPV3 

transient  receptor  potential 
cation  channel,  subfamily  V, 
member  3 

Plasma  Membrane 

ion  channel 

EFNA3 

ephrin-A3 

Plasma  Membrane 

kinase 

GUCY2C 

guanylate  cyclase  2C 

Plasma  Membrane 

kinase 

C16orf30 

chromosome  16  open  reading 
frame  30 

Plasma  Membrane 

other 

CCDC8 

coiled-coil  domain  containing  8 

Plasma  Membrane 

other 

EPB41L1 

erythrocyte  membrane  protein 
band  4.1-like  1 

Plasma  Membrane 

other 

VMD2L1 

bestrophin  2 

Plasma  Membrane 

other 

F7 

coagulation  factor  VII 

Plasma  Membrane 

peptidase 

NOTCH1 

Notch  homolog  1,  translocation- 
associated 

Plasma  Membrane 

transcription  regulator 

ASGR1 

asialoglycoprotein  receptor  1 

Plasma  Membrane 

transmembrane  receptor 

BRSK2 

BR  serine/threonine  kinase  2 

Unknown 

kinase 

C15orf39 

chromosome  15  open  reading 
frame  39 

Unknown 

other 

C19orf48 

chromosome  19  open  reading 
frame  48 

Unknown 

other 

C20orf166 

chromosome  20  open  reading 
frame  166 

Unknown 

other 

C21orf125 

chromosome  21  open  reading 
frame  125 

Unknown 

other 

CCDC124 

coiled-coil  domain  containing 

124 

Unknown 

other 

CRAMP1L 

Crm,  cramped-like 

Unknown 

other 

FBX044 

F-box  protein  44 

Unknown 

other 

KIAA1984 

KIAA1984 

Unknown 

other 

W81XWH-05- 1-0535;  PI:  Issa,  Annual  Report  2007,  Page  10 


KLHL30 

kelch-like  30 

Unknown 

other 

TBC1D16 

TBC1  domain  family,  member 

16 

Unknown 

other 

TMEM79 

transmembrane  protein  79 

Unknown 

other 

We  analyzed  the  functional  relationships  of  these  genes  by  Ingenuity  Pathway  Analysis 
software.  Two  major  networks  showed  significant  representation  of  the  genes  detected  by 
methylation  microarrays:  (1)  Cell  death,  cancer,  and  cell  proliferation  network  clustering 
around  Akt,  PI3K,  and  NFkB;  (Figure  1)  and  (2)  Cancer,  cell  death,  and  immunological 
disease  network  clustering  around  beta-estradiol,  progesterone,  and  TNF  (Figure  2). 
Twenty  seven  genes  from  the  top  50  list  could  be  classified  by  their  function.  Cell 
signaling  and  proliferation  were  the  most  represented  functions  (Table  2). 

Table  2.  Functional  classification  of  genes  methylated  in  MPD  patients. 


Function 

Genes 

Percent 

P  value 

Cell  Signaling 

11 

41% 

0.005 

Cellular  Growth  and  Proliferation 

10 

37% 

0.008 

Immunological  Disease 

6 

22% 

0.007 

Hematological  System  Development  and  Function 

7 

26% 

0.008 

Nucleic  Acid  Metabolism 

4 

15% 

0.006 

Cellular  Development 

11 

41% 

0.009 

Immune  and  Lymphatic  System  Development  and 

Function 

7 

26% 

0.008 

Developmental  Disorder 

5 

19% 

0.008 

Cell  Morphology 

7 

26% 

0.006 

Cell  Death 

6 

22% 

0.008 

Gene  Expression 

5 

19% 

0.007 

Cancer 

8 

30% 

0.009 

Cell  Cycle 

6 

22% 

0.009 

Cellular  Movement 

6 

22% 

0.008 

Cellular  Compromise 

3 

11% 

0.005 
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Figure  1.  Cell  death,  cancer,  and  cell  proliferation  network. 

Functional  relationships  were  prediceted  by  Ingenuity  Pathways  Analysis  software. 
Genes  detected  as  methylated  are  shaded.  Akt  kinase,  PI3  kinase,  and  NFkB  transcription 


factor  form  central  nodes. 
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Figure  2.  Cancer,  cell  death,  and  immunological  disease  network.  Functional 
relationships  were  prediceted  by  Ingenuity  Pathways  Analysis  software.  Genes  detected 
as  methylated  are  shaded.  MYC,  TNF,  beta-estradiol,  and  progesterone  form  central 


nodes. 
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Canonical  pathways  with  genes  hypermethylated  in  MPD  patients  are  shown  in  Table  3. 

Table  3.  Canonical  pathways  with  methylated  genes  in  MPD  patients. 


Canonical  Pathway 

Genes 

Purine  Metabolism 

MYH6,  GUCY2C,  PDE4C,  ADARB2 

Glycine,  Serine  and  Threonine  Metabolism 

SARDH,  CHDH 

Axonal  Guidance  Signaling 

EFNA3,  ABLIM2,  PRKCZ 

Cardiac  alpha-adrenergic  Signaling 

AKAP8,  PDE4C 

Insulin  Receptor  Signaling 

TSC2,  PRKCZ 

PI3K/AKT  Signaling 

TSC2,  PRKCZ 

cAMP-mediated  Signaling 

AKAP8,  PDE4C 

Synaptic  Long  Term  Depression 

GUCY2C,  PRKCZ 

NRF2-mediated  Oxidative  Stress  Response 

HSPB8,  PRKCZ 

HSPB8,  PRKCZ,  and  ME1  genes  detected  as  hypermethylated  participate  in  oxidative 
stress  pathways  (Table  4).  PRKCZ  gene  belongs  to  NFI-B  transcription  factor  pathway. 
NFIX  gene  is  a  member  of  the  same  family  as  NFI-B.  Interestingly,  NFI-B  upregulation 
was  reported  in  a  paper  showing  a  loss  of  heterozygosity  and  uniparental  disomy  in  the 
JAK2  locus  in  polycythemia  vera.15 

Table  4.  Metabolic  pathways  with  genes  hypermethylated  in  MPD  patients. 


Pathway 

Genes 

Oxidative  Stress  Response  Mediated  by  Nrf2  (NFE2L2) 

HSPB8,  PRKCZ 

Oxidative  Stress 

ME1 

Mechanism  of  Gene  Regulation  by  Peroxisome  Proliferators  via 

PPARa 

ME1 

NFI-B  Signaling  Pathway 

PRKCZ 

Fatty  Acid  Metabolism 

ME1 

Aryl  Hydrocarbon  Receptor  Signaling 

NFIX 
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Progress  on  Task  2 

Methylation  of  HOX  genes  in  MPD  patients. 

We  screened  31  HOX  genes  for  methylation  of  their  CpG  islands  in  PV  patients. 
Preliminary  results  were  presented  at  American  Association  for  Cancer  Research  Annual 
Meeting,  in  Los  Angeles,  CA,  in  April,  2007.  These  preliminary  results  suggested 
increased  methylation  in  HOXA4,  HOXC4,  HOXC9,  and  HOXC1 1  genes.  We  performed 
detailed  analysis  in  larger  numbers  of  MPD  patients  and  found  significantly  increased 
methylation  in  HOXA4,  HOXC4,  and  HOXC1 1  genes  (Figures  3-6). 
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Figure  3.  CpG  methylation  at  the  transcription  start  site  of  HOXA4  gene  in  blood 
cells  from  normal  controls  and  MPD  patients.  P<.0001,  Mann- Whitney  non- 


parametric  test. 
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HOXC4 


Figure  4.  CpG  methylation  at  the  transcription  start  site  of  HOXC4  gene  in  blood 
cells  from  normal  controls  and  MPD  patients.  P=.  01,  Mann- Whitney  non-parametric 
test. 
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Figure  5.  CpG  methylation  at  the  transcription  start  site  of  HOXC9  gene  in  blood 
cells  from  normal  controls  and  MPD  patients.  Although  the  difference  was  not 
statistically  significant,  about  25%  of  MPD  patients  showed  HOXC9  methylation  outside 


of  the  normal  range. 
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HOXC11 


Controls  (n=23)  MPD  (n=76) 


Figure  6.  CpG  methylation  at  the  transcription  start  site  of  HOXC11  gene  in  blood 
cells  from  normal  controls  and  MPD  patients.  P<.0001,  Mann- Whitney  non- 
parametric  test. 


Progress  on  Task  3 

Progesterone  receptor  functional  studies 

To  assess  the  functional  significance  of  progesterone  receptor  silencing,  we  explored  the 
effect  of  mifepristone,  a  PGR  antagonist,  on  in  vitro  response  of  BFU-E  erythroid 
progenitors  to  erythropoietin.  Mifepristone  treatment  of  normal  peripheral  blood 
increased  cultured  in  vitro  in  the  presence  of  erythropoietin  increased  the  expression  of 
genes  associated  with  erythroid  program:  transferrin  receptor,  glycophorin  A,  janus 


kinase  2,  and  beta  globin  (Fig.  7). 
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Figure  7.  Mifepristone  treatment  increases  the  expression  of  erythroid  genes. 

Normal  peripheral  blood  mononuclear  cells  were  cultured  in  vitro  for  14  days  in  the 
presence  of  erythropoietin  (3U/ml)  and  mifepristone  at  0,  10-7M  and  10-6M 
concentration.  Expression  of  mRNA  was  measured  by  real  time  PCR.  TFRC,  transferrin 
receptor,  EpoR,  erythropoietin  receptor,  GYPA,  glycophorin  A,  JAK2,  janus  2  kinase,  b- 
globin,  beta  globin  mRNA.  Fold  increase  of  mRNA  expression  over  cultures  without 
mifepristone  is  shown. 

Epigenetic  control  of  PRV-1  gene  expression 


We  found  that  expression  of  the  polycythemia  rubra  vera-1  (PRV-1)  gene  is  inversely 
correlated  with  DNA  methylation  status  of  CpG  sites  the  the  gene  transcription  start  site. 
PRV-1/CC177  is  a  GPI-linked  protein  that  is  expressed  on  a  subgroup  of  neutrophils. 
The  number  of  PRV-1  expressing  neutrophils  increases  in  pregnancy  and  sepsis,  or  after 
administration  of  granulocyte-colony  stimulating  factor.  Expression  of  the  PRV-1  gene 
is  also  increased  in  patients  with  polycythemia  vera  and  essential  thrombocythemia.  We 
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investigated  whether  DNA  methylation  of  the  PRV-1  gene  has  a  role  in  regulation  of 
transcription  and  expression  of  the  PRV-1  protein.  We  compared  the  level  of  methylation 
of  the  PRV-1  gene  and  expression  of  the  PRV-1  mRNA  in  nonnal  neutrophils  expressing 
PRV-1  to  those  that  are  PRV-1  negative.  We  also  studied  PRV-1  methylation  and 
mRNA  expression  in  patients  with  non-CML  myeloproliferative  disorders  and  in  an  in 
vitro  model  of  DNA  demethylation.  We  found  that  methylation  of  CpG  dinucleotides 
close  to  initiation  codon  of  the  PRV-1  gene  was  inversely  related  to  the  expression  of 
PRV-1  in  nonnal  neutrophils.  Furthermore,  overexpression  of  the  PRV-1  gene  in 
polycythemia  vera  (PV)  and  essential  thrombocythemia  (ET)  is  associated  with  a 
decrease  in  methylation  of  this  gene.  Among  patients  with  PV  and  ET,  methylation  of 
the  PRV-1  gene  is  also  inversely  correlated  with  the  presence  of  the  JAK2V617F  somatic 
mutation.  In  an  in  vitro  model,  exposure  of  KG1  and  KG  la  cells  to  a  DNA 
demethylating  agent  caused  a  decrease  in  methylation  of  the  PRV-1  gene  and  increased 
its  mRNA  level.  We  conclude  that  DNA  methylation  regulates  PRV-1  expression  under 
physiologic  and  pathologic  conditions.  The  results  were  published  in  Experimental 
Hematology  (Appendix  4) 

Effects  of  DNA  demethylating  treatment  in  MPD  patients. 

The  mechanism  of  DNA-demethylating  action  of  2-deoxy-5-aza-cytidine  (decitabine)  in 
vivo  is  not  fully  understood.  We  studied  the  dynamics  of  neoplastic  cell  clearance  during 
decitabine  treatment  (100  mg/m2/course  every  4  weeks)  using  quantitative  monitoring  of 
mutant  alleles  by  pyrosequencing  in  patients  with  chronic  myelomonocytic  leukemia 
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harboring  JAK2  or  NPM1  mutations.  CMML  patients  were  first  screened  for  JAK2  and 
NPM1  mutations,  and  three  patients  with  mutations  were  identified.  Mutant  allele 
percentages  in  mononuclear  cell  DNA  were  followed  after  treatment,  along  with 
methylation  of  LINE  1  and  ten  other  genes.  The  clearance  of  mutant  alleles  was  modest 
after  the  first  cycle,  despite  induction  of  hypomethylation.  Delayed  substantial  clearance 
was  observed  after  2-4  cycles  that  correlated  with  clinical  response.  Two  patients  had 
complete  disappearance  of  mutant  alleles  and  sustained  clinical  remissions.  In  another 
patient,  mutant  allele  was  detectable  at  clinical  remission,  which  lasted  for  8  months.  Our 
data  suggest  a  predominantly  non-cytotoxic  mechanism  of  action  for  decitabine,  leading 
to  altered  biology  of  the  neoplastic  clone  and/or  normal  cells.  The  result  are  in  press  in 
the  Blood  journal  (Appendix  5) 
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Key  Research  Accomplishments 


•  Used  DNA  oligonucleotide  microarrays  to  map  genome- wide  methylation  pattern 
in  MPD  patients 

•  Studied  DNA  methylation  dynamics  and  molecular  response  after  DNA- 
demethylation  therapy  in  MPD  patients. 

•  Characterized  epigenetic  control  by  DNA  methylation  of  the  polycythemia  rubra 
vera  gene  1  (PRV-1/CD177)  expression  on  neutrophils  of  MPD  patients  and 
normal  individuals. 

Reportable  Outcomes 

Meeting  presentations 

Poster,  DNA  methylation  of  HOX  genes  in  leukemia  and  myeloproliferative  disorders. 

American  Association  for  Cancer  Research  Annual  Meeting,  Los  Angeles,  CA,  April  14- 

18,2007. 

Meeting  abstracts  published 

Nussenzveig  RH,  Swierczek  S,  Jelinek  J,  Verstovsek  S,  Prchal  J,  Prchal  JT.  A  Novel  and 

Quantitative  Assay  To  Detect  JAK2  V617F  Allele  by  Real-Time  AS-PCR  and  Its 

Applicability  to  PV  Initiating  Mutation.  Blood  2006;  108:  3 1 5b-3 1 6b  (Appendix  1) 
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Oki  Y,  Jelinek  J,  Kantarjian  HM,  Issa  JPJ.  Hypomethylation  Induction  and  Molecular 
Response  after  Decitabine  Therapy  in  Chronic  Myelomonocytic  Leukemia  (CMML). 
Blood  2006;  108:  658a  (Appendix  2) 

Manuscript  published 

Nussenzveig  RH,  Swierczek  SI,  Jelinek  J,  Gaikwad  A,  Liu  E,  Verstovsek  S,  Prchal  JF, 
Prchal  JT.  Polycythemia  vera  is  not  initiated  by  JAK2V617F  mutation.  Experimental 
Hematology.  2007;35:32-38.  (Appendix  3) 

Jelinek  J,  Fi  J,  Mnjoyan  Z,  Issa  JP,  Prchal  JT,  Afshar-Kharghan  V.  Epigenetic  control  of 
PRV-1  expression  on  neutrophils.  Exp  Hematol.  2007;35(1 1):  1677-1683.  (Appendix  4) 

Manuscript  in  press 

Yasuhiro  Oki*,  Jaroslav  Jelinek*,  Lanlan  Shen,  Hagop  M.  Kantarjian,  and  Jean-Pierre  J. 
Issa.  Induction  of  hypomethylation  and  molecular  response  after  decitabine  therapy  in 
patients  with  chronic  myelomonocytic  leukemia.  Blood,  in  press.  (Appendix  5) 

Manuscript  in  preparation 

Jelinek  J,  He  R,  Bueso-Ramos  CE,  Verstovsek  S,  Prchal  JT,  Issa  JPJ.  Methylation  of 
progesterone  receptor  CpG  islands  in  polycythemia  vera,  myelofibrosis  and  leukemia. 


W81XWH-05- 1-0535;  PI:  Issa,  Annual  Report  2007,  Page  22 


CONCLUSIONS 


Our  data  show  that  PGR  and  CDH13  CpG  islands  are  methylated  in  a  significant 
proportion  of  PV  and  MF  patients  and  in  leukemia.  Silencing  of  these  genes  by 
methylation  may  contribute  to  disease  development  by  altering  the  response  of 
hematopoietic  cells  to  proliferative  stimuli  or  their  interactions  with  stroma. 


“So  what:” 


Epigenetic  silencing  by  cytosine  methylation  in  selective  CpG  islands  may  play  a  role  in 
the  development  of  myeloproliferative  disorders.  The  hypomethylating  drug  decitabine 
may  be  considered  for  clinical  trials  in  patients  non-responding  to  conventional  treatment. 
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Appendix  1 


A  Novel  and  Quantitative  Assay  To  Detect  JAK2  V617F 
Allele  by  Real-Time  AS-PCR  and  Its  Applicability  to 
PV  Initiating  Mutation. 

Roberto  H.  Nussenzveig,  PhD1’*,  Sabina  Swierczek,  PhD2’*,  Jaroslav  Jelinek,  MD1, 
Srdan  Verstovsek,  MD1,  Jaroslav  Prchal,  MD3  and  Josef  T.  Prchal,  MD2 

1  Dept.  Leukemia,  MD  Anderson  Cancer  Center,  Houston,  TX,  USA; 2  Division  of  Hematology,  SOM, 
University  of  Utah,  Salt  Lake  City,  UT,  USA  and  3  Dept.  Oncology,  McGill  University,  Montreal,  QC, 
Canada. 


Abstract 

Polycythemia  vera  (PV)  arises  due  to  a  somatic  mutation(s)  of  a  single  hematopoietic 
stem  cell  leading  to  clonal  hematopoiesis.  Greater  than  80%  of  PV  patients  carry  a 
somatic  mutation  in  JAK2  (V617F).  Growing  evidence  suggests  that  increased  frequency 
of  the  JAK2V617F  allele  may  have  a  prognostic  impact  on  certain  clinical  aspects  of  PV, 
and,  possibly,  in  other  myeloproliferative  disorders  associated  with  this  mutation.  We 
have  developed  a  novel  approach  to  primer  design  for  Real-Time  quantitative  allele- 
specific  PCR.  Allelic  discrimination  is  enhanced  by  the  combined  synergistic  effects  of  an 
artificial  mismatch  introduced  in  the  -1  position,  starting  from  the  3’  end  of  the  primer, 
and  the  use  of  a  locked  nucleic  acid  (LNA)  modified  nucleoside  placed  at  the  -2  position. 
We  provide  evidence  that  the  -2  LNA  assists  in  stabilizing  the  3’  end,  while  the  -1 
mismatch  provides  specificity  but  not  stability.  The  difference  in  cycle  number  between 
the  two  allele-specific  reactions  is  used  to  calculate  the  relative  allele  frequencies.  We 
demonstrate  the  robustness,  sensitivity  and  reproducibility  of  our  design.  The  proportion 
of  mutant  JAK2  allele  detennined  by  pyrosequencing  and  kinetic  allele-specific  PCR  was 
highly  concordant  with  an  average  allele  frequency  deviation  of  2.6%.  Repeated 
determination  of  allelic  ratios  in  multiple  patient  samples  was  highly  reproducible  with  a 
standard  deviation  of  1.5%.  We  have  also  determined  that  the  design  and  assay  is  highly 
sensitive;  as  little  as  0. 1%  mutant  allele  in  40-50  ng  of  genomic  DNA  can  be  detected. 

We  further  tested  the  applicability  of  this  technique  to  the  analysis  of  individual  BFU-E 
colonies  in  order  to  address  the  question  whether  the  JAK2V617F  is  the  disease  initiating 
mutation.  Less  than  10%  of  a  single  isolated  BFU-E  colony,  originating  from  a  single 
progenitor,  is  sufficient  for  determination  of  allele  frequency.  The  remainder  of  the 
colony  may  be  used  for  other  analyses.  A  proportion  of  0  or  50  or  100  percent  JAK2 
mutant  allele  is  expected  from  each  individual  BFU-E  colony,  which  was  indeed 
observed.  However,  when  we  tested  granulocytes  from  PV  females,  wherein  the 
granulocytes  were  found  to  be  clonal  by  the  X-chromosome  transcriptionally  based 
clonality  assay,  we  found  3  females  <50  (27.5  ±1 1)  and  7  females  >50  (75  ±10.5)percent 
mutant  JAK2  allele  frequencies.  This  result  suggests  the  presence  of  a  heterogeneous 
population  of  cells  with  differing  genotypes  regarding  the  JAK2  mutant  allele,  and  is 
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further  supported  by  our  genotyping  results  with  individual  BFU-E  colonies  as  described 
above.  Our  PV  data  suggest  that  the  JAK2V617F  may  not  be  the  PV  initiating  mutation. 
This  novel  primer  design  is  simple,  does  not  require  tedious  optimization  of  reaction 
conditions,  and  can  be  applied  to  any  kinetic  PCR platform  for  reliable  and  sensitive 
determination  of  allele  frequencies.  Potential  applications  are  varied,  such  as,  quantitative 
determination  of  mosaicism,  proportion  of  fetal  cells  in  maternal  circulation,  detection  of 
minimal  residual  disease  associated  with  known  somatic  mutation  (such  as  reduction  of 
malignant  cells  by  chemotherapy  or  reappearance  of  resistant  clone),  rapid  monitoring  of 
efficacy  of  new  drugs  in  both  "in  vitro"  systems  as  well  as  clinical  trials,  and  many  others 
that  require  quantitation  of  allele  frequencies. 

Appendix  2 


Hypomethylation  Induction  and  Molecular  Response 
after  Decitabine  Therapy  in  Chronic  Myelomonocytic 
Leukemia  (CMML). 

Yasuhiro  Oki,  MD*’1’2,  Jaroslav  Jelinek,  M.D.  Ph.D.1,  Hagop  M.  Kantarjian,  M.D.1 
and  Jean-Pierre  J.  Issa,  M.D.1 

1  Department  of  Leukemia,  University  of  Texas  M.D.  Anderson  Cancer  Center,  Houston,  TX,  USA  and  2 
Department  of  Clinical  Oncology  and  Hematology,  Aichi  Cancer  Center,  Nagoya,  Aichi,  Japan. 

Abstract 

Decitabine  has  shown  therapeutic  activity  in  patients  with  MDS  and  CMML.  The 
mechanisms  of  response  to  therapy  remain  incompletely  understood.  In  particular,  the 
relative  contribution  of  this  drug’s  ability  to  induce  hypomethylation  and  cytotoxicity 
remains  unclear.  To  address  this  issue,  we  studied  the  dynamics  of  neoplastic  cell 
clearance  during  decitabine  treatment  determined  by  quantitative  monitoring  of  the 
mutant  allele  using  pyrosequencing.  DNA  extracted  from  peripheral  blood  mononuclear 
cells  from  consented  patients  with  CMML  in  a  decitabine  phase  II  study  were  first 
screened  for  JAK2  and  NPM1  mutations  as  previously  reported.  We  identified  three 
patients  with  mutations  (two  with  JAK2  mutation,  one  with  NPM 1  mutation)  and  samples 
at  multiple  points  during  therapy  were  available.  All  three  carried  nonnal  karyotype. 

LINE  repetitive  element  methylation  and  several  other  gene  specific  methylations  were 
also  assessed.  In  the  three  patients,  LINE  methylation  decreased  after  each  cycle  of 
therapy,  and  recovered  to  near  baseline  after  the  drug  was  stopped  (e.g.  during  the  first 
cycle,  average  relative  hypomethylation  from  baseline  was  13.9%  at  day  12  and  6.5%  at 
day  28).  At  the  same  time,  the  proportion  of  circulating  neoplastic  cells  decreased  slowly 
after  the  first  cycle  (decrease  by  19.3%  at  day  12  and  13.5%  at  day  28).  A  substantial 
decrease  in  mutant  allele  percentage  was  observed  after  cycles  2,  3,  and  2  in  patients  1,  2, 
and  3,  respectively.  Clinical  complete  responses  were  achieved  along  with  molecular 
responses  at  cycles  5,  4  and  2,  respectively.  Patients  1  and  2  showed  complete 
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disappearance  of  detectable  neoplastic  clones,  and  had  sustained  remissions  (duration  1.5 
and  2.5  years).  In  patient  3,  the  proportion  of  neoplastic  cells  was  lower  than  baseline  but 
still  detectable  at  clinical  remission,  and  the  remission  only  lasted  8  months.  We  conclude 
that  neoplastic  cell  clearance  after  decitabine  therapy  in  CMML  is  observed  after  several 
courses  of  therapy,  and  is  initially  seen  concurrently  with  hypomethylation.  While  LINE 
methylation  returns  to  its  steady  state  values  after  completion  of  decitabine  infusion,  the 
tumor  elimination  process  slowly  continues.  Our  data  suggest  a  non-cyto toxic  mechanism 
of  action  for  the  drug,  whereby  the  biology  of  the  neoplastic  clone  is  altered  by 
hypomethylation,  leading  to  delayed  clearances  of  unknown  mechanism.  Possibilities 
include  an  immune  response  and  effects  on  the  neoplastic  (or  normal)  stem  cells. 


Page  28 


ELSEVIER 


EXPERIMENTAL 

HEMATOLOGY 


Experimental  Hematology  35  (2007)  32-38 


Polycythemia  vera  is  not  initiated  by  JAK2V617F  mutation 
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Objective.  The  somatic  JAK2V6,7F  mutation  is  seen  in  most  polycythemia  vera  (PV)  patients; 
however,  it  is  not  clear  if  JAK2V6l7F  is  the  PV-initiating  mutation. 

Methods.  In  order  to  examine  this  issue,  we  developed  a  novel  real-time  quantitative  allele- 
specific  PCR,  in  which  allelic  discrimination  is  enhanced  by  the  synergistic  effect  of  a  mis¬ 
match  in  the  —1  position,  and  a  locked  nucleic  acid  (LNA)  nucleoside  at  the  —2  position. 
Results.  Determination  of  allelic  frequencies  was  reproducible  (SD  =  1.5%)  and  sensi¬ 
tive— 0.1%  mutant  allele  detected  in  40  ng  of  DNA.  The  JAK2V617F  frequency  in  clonal  gran¬ 
ulocytes  from  3  PV  females  was  less  than  50%  (27.5  ±11)  and  in  7  females  greater  than  50% 
(75  ±  10.5).  We  also  found  that  wild-type  JAK2  BFU-E  colonies  from  PV  patients  can  grow 
without  erythropoietin.  The  identification  of  the  primary  genetic  lesion  resulting  in  PV  is 
essential  for  the  development  of  novel  therapeutic  strategies. 

Conclusion.  Our  studies  correlating  the  frequency  of  JAK2V6llv  mutant  allele  and  clonality, 
as  well  as  the  presence  of  homozygous  wild-type  JAK2  erythropoietin-independent  erythroid 
colonies,  provide  compelling  evidence  that  the  JAK2V6m  is  not  the  PV-initiating  mutation. 
This  supports  a  model  wherein  the  JAK2yt’l7F  mutation  arises  as  a  secondary  genetic  event. 
Furthermore,  our  results  indicate  that  an  undefined  molecular  lesion,  preceding  JAK2V6i7r, 
is  responsible  for  clonal  hematopoiesis  in  PV.  We  conclude  that  development  of  therapeutic 
strategies  that  target  the  JAK2V6I7F  clonal  cells  may  not  be  sufficient  for  eradication  of 
PV.  ©  2007  International  Society  for  Experimental  Hematology.  Published  by  Elsevier  Inc. 


Myeloproliferative  disorders  (MPDs)  are  hematological 
malignancies  due  to  clonal  proliferation  arising  from  a  sin¬ 
gle  multilineage  stem  cell.  The  World  Health  Organization 
(WHO)  classified  MPDs  into  the  following  four  prototypi¬ 
cal  clinical  diseases:  polycythemia  vera  (PV),  essential 
thrombocythemia  (ET),  myelofibrosis  with  myeloid  meta¬ 
plasia  (MMM),  and  chronic  myeloid  leukemia  (CML).  Sub¬ 
sequently,  other  related  disorders  were  added  to  this 
classification  [1],  The  chromosomal  translocation  t(9,22), 
resulting  in  a  chimeric  fusion  protein  encoded  by  the 
Bcr-Abl  gene,  leads  to  proliferation  and  survival  of  myeloid 
progenitor  cells  in  CML.  Recently,  a  novel  somatic  single¬ 
point  mutation  in  the  9p  chromosomal  region  encoding  the 
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tyrosine  kinase  JAK2  (1849  G  to  T)  has  been  reported  in 
several  MPDs  [2-6].  This  mutation  results  in  a  nonsynony- 
mous  amino  acid  substitution  at  position  617  (valine  to 
phenylalanine)  located  in  the  JH2  pseudo-kinase  auto- 
inhibitory  domain  [7],  The  mutation  renders  the  enzyme 
constitutively  active  and  leads  to  cytokine  hypersensitivity 
and  erythrocytosis  in  a  mouse  model  [8].  The  JAK2V6nf 
mutation  has  been  reported  in  greater  than  80%  of  PV  pa¬ 
tients  and  in  approximately  half  of  ET  or  MMM  patients 
[5].  The  proportion  of  the  mutant  allele  is  highly  variable, 
with  approximately  30%  of  PV  patients  having  loss-of-het- 
erozygosity  of  chromosome  9p  created  in  most  instances  by 
uniparental  disomy  [9], 

Studies  in  MPD  support  that  acquisition  of  the 
JAK2x'6,7F  mutation  results  from  a  somatic  event.  The  pro¬ 
portion  of  JAK2  mutant  allele  in  a  clonal  granulocyte  cell 
population  from  the  peripheral  blood  of  a  patient  can  vary 
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from  0  to  100%.  According  to  one  scenario,  the  mutation  is 
not  the  first  event  leading  to  disease;  this  is  further  sup¬ 
ported  by  the  existence  of  PV  families  with  a  predisposition 
to  acquire  PV  [10].  Moreover,  Kralovics  and  colleagues 
[11]  report  the  identification  of  a  PV  patient  with  1% 
JAK2  mutant  allele  frequency,  yet  harboring  the  character¬ 
istic  MPD  del20q  in  100%  of  her  granulocytes.  To  date,  in 
spite  of  the  identification  of  a  single  patient,  this  report  pro¬ 
vides  the  most  compelling  evidence  that  JAK2V6I7F  is  not 
the  primary  event  leading  to  PV.  In  addition,  zygosity  in  pa¬ 
tients  carrying  the  JAK2V6I7F  mutation  can  be  heteroge¬ 
neous,  with  wild-type,  heterozygous,  and  mutant 
homozygous  clones  being  present.  The  clinical  significance 
and  role  of  the  JAK2  mutation  in  the  pathology  and  progres¬ 
sion  of  disease,  especially  for  PV,  is  being  investigated.  It 
appears  that  some  complications  of  MPDs  such  as  the  de¬ 
gree  of  marrow  fibrosis  or  thrombotic  tendencies  correlate 
with  overall  proportion  of  the  mutant  allele  in  circulating 
cells  [8,12].  The  close  association  of  JAK2WM1F  with  PV 
and  burst-forming  unit  erythroid  (BFU-E)  endogenous  col¬ 
onies  of  PV  makes  this  mutation  an  excellent  disease 
marker.  Hence,  disease  progression,  treatment,  and  devel¬ 
opment  of  novel  chemotherapeutics  can  be  readily  assessed 
by  following  changes  in  the  proportion  of  the  mutant  JAK2 
allele  in  clonal  granulocytes. 

To  date,  several  different  methods  have  been  developed 
to  quantitatively  determine  the  proportion  of  single  nucleo¬ 
tide  polymorphisms  (SNPs)  in  pooled  DNA  samples  [13]. 
These  include  (but  are  not  limited  to)  allele-specific  PCR 
(AS-PCR)  and  pyrosequencing.  A  number  of  strategies 
have  been  employed  to  improve  the  specificity  and  reliabil¬ 
ity  of  AS-PCR  and  adapt  it  to  real-time  monitoring.  These 
modifications  have  demonstrated  that  AS-PCR  can  be  a  re¬ 
liable  tool  for  genotyping  provided  time  is  taken  to  care¬ 
fully  design  and  optimize  PCR  conditions.  On  the  other 
hand,  pyrosequencing  requires  instrumentation  that  is  not 
readily  available  and  setup  costs  are  prohibitive  for  routine 
laboratory  testing.  Moreover,  pyrosequencing  cannot  be 
used  for  routine  screening  of  residual  or  early  disease,  since 
the  proportion  of  the  JAK2  mutant  allele  is  below  the  5% 
detection  limit  threshold. 

We  describe  here  the  development  of  a  novel  approach 
to  quantitatively  determine  the  proportion  of  wild-type 
(wt)  and  mutant  JAK2  alleles  by  a  novel  modification  of 
real-time  AS-PCR.  We  demonstrate  the  high  sensitivity, 
specificity,  and  reproducibility  of  our  method  that  permits 
accurate  quantitation  of  the  JAK2V6I7F  mutant  allele  using 
now  widely  available  kinetic  PCR  instrumentation.  More¬ 
over,  our  approach  to  allele-specific  primer  design  is  sim¬ 
ple,  does  not  require  optimization  of  PCR  conditions,  and 
incorporates  nucleotides  that  are  readily  available  from 
most  oligonucleotide  synthesis  companies.  We  used  this 
methodology  to  enhance  our  understanding  of  the  molecu¬ 
lar  biology  of  PV.  Here,  we  provide  definitive  evidence  that 
the  JAK2V6I7F  mutation  is  not  the  PV-initiating  event. 


Materials  and  methods 

Samples 

Blood  samples  were  collected  from  20  unrelated  PV  patients  and  3 
healthy  volunteers  using  an  Institutional  Review  Board-approved 
protocol.  Granulocytes  were  isolated  from  peripheral  blood  using 
standard  methods.  Genomic  DNA  (gDNA)  was  extracted  from  pe¬ 
ripheral  blood  leukocytes  and  purified  granulocytes  using  the 
Puregene  DNA  purification  kit  as  recommended  by  the  manufac¬ 
turer  (Gentra,  Minneapolis,  MN,  USA). 

BFU-E  colony  culture 

In  vitro  assay  of  erythroid  progenitors’  responsiveness  to  erythro¬ 
poietin  (EPO)  was  performed  as  previously  described  [14]. 
Briefly,  mononuclear  cells  from  peripheral  blood  were  isolated 
on  Histopaque  1077  (Sigma,  St.  Louis,  MO,  USA)  and  cultured 
at  a  final  density  of  3  x  10s  cells/mL  in  Methocult  H-4531  me¬ 
dium  (StemCell  Technologies  Inc.,  Vancouver,  BC,  Canada)  in 
35-mm  Petri  dishes  in  the  presence  or  absence  of  3  U/mL  EPO. 
Cultures  were  maintained  at  37°C  and  5%  carbon  dioxide.  Single 
erythroid  colonies,  selected  using  standard  criteria,  were  picked 
after  14  days  in  culture  using  micropipettes. 

Real-time  AS-PCR 

Real-time  AS-PCR  was  performed  on  an  Applied  Biosystems 
7000  Sequence  Detection  System  (Applied  Biosystems,  Foster 
City,  CA,  USA).  Typical  reactions  (25  pL)  consisted  of  1  x  Taq- 
Man  Universal  PCR  master  mix  (Applied  Biosystems,  Foster 
City,  CA,  USA);  300  nM  JAK2  universal  forward  and  allele-spe¬ 
cific  reverse  primers;  125  nM  FAM-labeled  JAK2  MGBNFQ 
probe  (Applied  Biosystems);  and  1  to  50  ng  purified  genomic 
DNA.  Enzyme  activation  (95°C  for  10  minutes)  was  followed 
by  50  cycles  of  95°C  for  15  seconds,  and  60°C  for  I  minute.  Al¬ 
lele-specific  primers  were  designed  using  the  software  program 
Oligo  6.7  (Molecular  Biology  Insights,  Inc.,  Cascade,  CO,  USA) 
and  synthesized  by  IDT  (Integrated  DNA  Technologies,  Coral- 
ville,  IA,  USA).  The  3'  terminal  sequence  of  the  reverse  primer 
(C  or  A)-gJAK2-R  was  selected  for  specific  amplification  of  the 
mutant  allele.  Allelic  discrimination  was  enhanced  by  introducing 
an  artificial  mismatch  (T:T)  in  the  -  I  position  starting  from  the  3' 
end  of  the  primer,  as  described  in  the  amplification  refractory  mu¬ 
tation  system  (ARMS,  [15]).  An  additional  modified  locked  nu¬ 
cleic  acid  (LNA)  base  [16]  was  placed  at  the  -2  position  (G). 
Control  primers  either  with  modification  alone  or  in  the  absence 
of  any  modification  were  tested  as  described  above.  Sequences 
and  specifications  of  ail  printers  and  the  detection  probe  are 
provided  as  supplemental  Table  S 1 . 

Pyrosequencing 

Quantitation  of  the  JAK2V617F  mutant  allele  was  done  as  previ¬ 
ously  described  [17]  using  a  PSQ  HS  96  pyrosequencer  and  re¬ 
agents  provided  by  the  manufacturer  (Biotage,  Uppsala,  Sweden). 

Allele  frequency  calculations 

Allele  frequency  was  calculated  as  described  by  Germer  et  al. 
[13],  The  difference  in  cycle  threshold  (AC,)  between  the  two  al¬ 
lele-specific  PCR  reactions  is  a  measure  of  the  proportion  or  fre¬ 
quency  of  the  allele,  assuming  that  initial  replication  efficiency  is 
100%.  If  the  amplification  efficiencies  of  the  two  allele-specific 
reactions  differ  slightly,  this  can  be  corrected  by  measuring  the 
AC,  on  a  DNA  sample  known  to  be  heterozygous  for  the  mutation 
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Figure  1.  Comparison  of  primer  modifications  for  discrimination  of  wild-type  and  mutant  JAK2V6‘71'  alleles.  Genomic  DNA  from  a  healthy  donor  (JAK2 
homozygous  wild-type)  and  a  PV  patient  ( JAK2  homozygous  mutant)  were  used  to  evaluate  the  performance  of  G-  and  T-allele-specific  primers  bearing 
different  permutations  of  mismatch  and  LNA  modifications,  and  with  no  modifications  (see  supplemental  Table  SI).  (A)  ACt  for  G-allele-specific  primers; 
more  negative  indicates  greater  discrimination.  (B)  ACt  for  T-allele-specific  primers;  more  negative  indicates  greater  discrimination.  (C)  Ct  delay  resulting 
from  introduction  of  modification  in  G-allele-specific  primers  compared  to  G-primer  without  modification.  (D)  Ct  delay  resulting  from  introduction  of  mod¬ 
ification  in  T-allele-specific  primers  compared  to  T-primer  without  modification.  Data  represent  mean  ±  SD  of  three  independent  determinations. 


of  interest.  The  AC,  in  the  heterozygous  sample  should  be  zero; 
any  deviation  from  zero  can  be  subtracted  from  all  AC,  measure¬ 
ments  in  order  to  compensate  for  differing  amplification  efficien¬ 
cies  and  is  represented  by  IIC  (heterozygote  correction  factor). 
Therefore,  HCAC,  represents  the  heterozygote  corrected  differ¬ 
ence  in  cycle  threshold  between  the  two  allele-specific  PCR 
reactions. 

AC,  =  C, ^allele,  —  C,-allele2  (1) 

and 

HCAC,  =  AC,  -  (HC  C, .allele,  -  HC  C,.allele2)  (2) 

Results  obtained  in  equations  1  and  2  are  used  to  find  the  fre¬ 
quency  of  allele;  in  equation  3: 


Freq. allele,  =  1  /  ( EHCACl  +1)  (3) 

where  E  represents  the  efficiency  of  PCR  amplification  for 
allelei  and  can  be  deduced  by  the  slope  of  serially  diluted 
sample. 


Statistical  calculations 

One-way  ANOVA  and  multiple  comparison  (Newman-Keuls  test) 
were  used  to  analyze  statistically  significant  differences  between 
ACt  and  delay  results  (Fig.  1)  when  comparing  primer  design 
modifications.  Statistical  significance  was  assumed  for  p  <  0.05. 
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Results 

Assay  design 

Allele-specific  PCR  is  widely  used  for  SNP  genotyping  and 
is  based  on  amplification  of  DNA  by  an  allele-specific 
primer  matching  the  polymorphism  at  the  3'  position.  In 
theory,  the  allele-specific  primer  containing  the  mismatched 
nucleoside  at  the  3'  end  should  not  be  extended  by  Taq 
DNA  polymerase  (Taq).  However,  it  has  been  shown  that 
Taq  can  extend  mismatched  allele-specific  primers,  gener¬ 
ating  false-positive  results.  A  number  of  different  strategies 
have  been  developed  to  improve  specificity  and  reliability 
of  this  technique  [15,18-21]. 

Our  assay  has  made  use  of  two  such  strategies,  the  inclu¬ 
sion  of  a  second  mismatch  and  a  modified  LNA  base.  We 
hypothesized  that  inclusion  of  a  second  mismatch  at  the 
-1  position  would  further  destabilize  the  3'  end  of  the  al¬ 
lele-specific  mismatched  primer,  resulting  in  an  increase 
in  specificity  of  the  matched  primer  after  the  first  cycle.  Ol¬ 
igonucleotides  carrying  LNA  bases  are  thought  to  favor  the 
formation  of  A-helix  DNA  duplexes,  improved  base-pair 
stacking,  and  higher  melting  temperature  (Tm)  [16].  There¬ 
fore,  we  hypothesized  that  the  LNA  at  the  -2  position 
would  stabilize  the  matched  bases  of  our  allele-specific 
primer,  increase  primer  Tm,  enhance  allele-specific  ampli¬ 
fication,  and  decrease  assay-to-assay  variability.  A  sche¬ 
matic  representation  of  our  assay  design  is  provided  as 
Supplemental  Figure  SI. 

Real-time  AS-PCR 

Genomic  DNA  was  isolated  from  granulocytes  of  a  control 
(Ctl)  and  one  PV  patient  (PV4)  representing  the  two  homo¬ 
zygous  JAK2  genotypes  (GG  and  TT)  as  previously  deter¬ 
mined  by  pyrosequencing.  Comparison  of  our  novel 
allele-specific  primer  design,  with  control  allele-specific 
primers  bearing  a  mismatch,  LNA  modified  base,  and  no 
modification  (supplemental  Table  SI),  was  evaluated  by 
real-time  PCR  using  both  homozygous  genotypes  (Fig.  1). 
A  dynamic  increase  in  fluorescence,  reflecting  probe  cleav¬ 
age,  was  seen  with  all  primer  designs.  In  the  absence  of  any 
modification,  the  G-aliele-specific  primer  efficiently  ex¬ 
tended  both  matched  and  mismatched  targets,  resulting  in 
poor  discrimination  between  genotypes  (Fig.  1A).  Allele- 
specific  primers  carrying  the  LNA  modified  base  showed 
only  marginal  improvement  discriminating  between  geno¬ 
types  (Fig.  1  A,B).  Results  for  the  primers  carrying  the  extra 
mismatch  and  our  design  (Fig.  1A,B)  were  excellent  for 
both  wild-type  and  mutant  homozygous  genotypes,  with 
ACt  greater  than  14  cycles  between  matched  and  mis¬ 
matched  primers.  This  indicated  that  the  additional  mis¬ 
match,  in  this  particular  application,  was  responsible  for 
improved  discrimination  and  decreased  mismatch  extension 
by  the  G-allele-specific  primer. 

Next  we  examined  the  reproducibility  of  the  results  by 
determining  the  ACt  of  matched  and  mismatched  primers 


in  three  independent  measurements  (Fig.  1A,B).  Results 
indicate  that,  although  discrimination  is  significantly 
improved  by  the  introduction  of  the  mismatch  at  the  -1 
position  and  comparable  to  our  primer  design,  addition  of 
the  LNA  base  at  the  —2  position  enhances  reproducibility 
as  expected.  We  also  found  that  positioning  of  the  LNA 
base  at  either  the  —  1  or  —4  positions  slightly  improved  dis¬ 
crimination,  albeit  at  a  cost  of  decreasing  the  sensitivity  of 
the  assay  (Fig.  1C,D).  ACt  differences  obtained  for  the  G- 
ailele-specific  primers,  bearing  either  the  mismatch  alone 
(gJAK2-(-lmut)-R)  or  mismatch  and  —2  LNA  (gJAK2-(- 
2)LNA-R)  modifications,  were  not  significant  (p  >  0.05). 
In  contrast,  the  differences  in  ACt  for  the  T-allele-specific 
primers,  harboring  the  same  modifications  as  above,  were 
significant  (p  <  0.05). 

Serial  dilutions  of  genomic  DNA  from  both  wild-type 
and  mutant  homozygous  samples  were  assayed  by  real¬ 
time  PCR  using  our  allele-specific  primers.  PCR  reaction 
efficiency  was  determined  for  both  allele-specific  primers 
and  was  found  to  be  greater  than  92%.  Allele  detection 
limits  were  estimated  to  be  in  the  range  of  50  pg  genomic 
DNA  (approximately  7  diploid  genome  equivalents).  Since 
JAK2V6>7F  is  a  somatic  mutation,  we  determined  the  un¬ 
equal  amplification  efficiency  and  calculated  the  heterozy¬ 
gosity  correction  factor  (HC)  between  our  allele-specific 
primers  by  measuring  the  ACt  of  a  50:50  mix  of  our 
wild-type  and  homozygous  mutant  DNA  samples. 

Real-time  vs  pyrosequencing 

Since  pyrosequencing  is  considered  a  quantitative  and  reli¬ 
able  technique  for  the  detection  of  SNP,  we  decided  to  per¬ 
form  a  blind  comparison  between  our  allele  frequency 
determinations  and  those  independently  obtained  by  pyro¬ 
sequencing.  Table  1  shows  the  estimated  frequencies  of 
the  mutant  T-allete  in  granulocytes  isolated  from  different 
PV  patients  obtained  by  both  real-time  PCR  and  pyrose¬ 
quencing  with  an  excellent  correlation  in  a  linear  regression 
model  (y  =  1.01*  +  2.6;  R2  =  0.99;  p  <  0.0001)  (Supple¬ 
mental  Fig.  S2). 

Reproducibility  of  the  assay 

In  order  to  test  the  robustness  of  our  design,  the  frequency 
of  the  JAK2  mutant  T-allele  was  determined  in  31  genomic 
DNA  samples  from  peripheral  blood  and  granulocytes  of 
PV  patients  (Table  2).  Three  independent  measurements 
were  performed  for  each  sample.  Mean  allele  frequency 
was  found  to  deviate  by  less  than  1.5%,  on  average,  for 
all  samples  tested  (Table  2). 

Mutant  allele  frequency  sensitivity 
The  JAK2  mutant  allele  detection  sensitivity  was  investi¬ 
gated  by  mixing  normal  control  genomic  DNA  with  de¬ 
creasing  proportions  of  genomic  DNA  isolated  from  the 
human  erythro-leukemic  cell  line  (HEL),  which  is  homozy¬ 
gous  for  the  JAK2V617F  mutation  (Fig.  2).  The  most  sensi¬ 
tive  techniques  used  to  detect  low  levels  of  SNP  are  still 
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Table  1.  Comparison  of  T-allelc  frequency  obtained  by  pyrosequencing 
and  our  method 


Patient 

number 

Date  on  sample 

T-allele  (%) 
(pyrosequencing) 

T-allele  (%) 
(AS-PCR) 

PV1 

1 1/3/2003 

57.0 

64.8 

PV2 

3/1/2005 

78.3 

80.7 

PV3 

4/13/2005 

21.3 

23.4 

PV3 

6/30/2005 

28.7 

30.0 

PV4 

12/29/2004 

69.4 

66.4 

PV4 

3/1/2005 

59.2 

60.4 

PV4 

5/26/2005 

57.9 

68.4 

PV4 

6/30/2005 

18.8 

24.1 

PV5 

8/2/2005 

44.9 

46.6 

PV6 

8/31/2005 

0.0 

0.0 

PV7 

1/4/2005 

84.1 

88.5 

PV7 

3/30/2005 

66.7 

73.4 

PV7 

6/30/2005 

80.2 

82.2 

PV7 

7/12/2005 

78.4 

84.1 

PV7 

10/14/2005 

79.0 

81.4 

PV8 

10/19/2005 

96.0 

97.0 

PV9 

8/22/2005 

83.0 

87.3 

PV10 

5/11/2005 

80.0 

84.2 

limited  by  the  mass  of  gDNA  assayed.  The  mass  of  a  single 
haploid  human  genome  is  estimated  to  be  3.7  pg.  Therefore, 
if  40  ng  of  genomic  DNA  are  used  to  assay  for  the  presence 
of  a  mutated  allele  with  a  frequency  of  0.01%,  only  a  single 
copy  may  be  detected  by  random  chance.  Most  assays  use 
20  to  50  ng  gDNA;  hence,  a  reasonable  detection  limit 
should  be  0.1%.  We  have  routinely  and  reliably  measured 
JAK2  mutant  allele  at  a  frequency  of  0.1%  using  40  to  50 
ng  gDNA  (Fig.  2).  However,  even  though  sensitivity  of 
our  assay  has  not  been  determined  with  0.5  pg  of  gDNA, 
we  should  be  able  to  reliably  detect  less  than  0.01%  mutant 
allele  (Fig.  2). 

Is  JAK2V6l7h  the  disease-initiating  mutation? 

We  have  demonstrated  the  relevance  and  usefulness  of  our 
primer  design  while  investigating  mutant  JAK2  T-allele  fre¬ 
quency  and  clonality  in  several  female  PV  patients.  Geno¬ 
mic  DNA  isolated  from  granulocytes  of  10  female  PV 
patients,  who  were  found  to  be  clonal  by  the  X-chromo- 
some  transcriptional  clonality  assay  [22],  were  used  for 
the  determination  of  JAK2  mutant  allele  frequency.  In  3 
of  the  10  female  patients  with  clonal  granulopoiesis,  mutant 
allele  frequencies  were  below  50%  (27.5  ±11)  while  the 
remaining  7  had  frequencies  greater  than  50%  (75  ± 
10.5).  This  result,  in  combination  with  the  clonality  data, 
is  indicative  of  a  heterogeneous  population  composed  of 
granulocytes  with  differing  JAK2V6I7F  genotypes  (GG, 
GT,  or  TT)  and  reinforces  the  model  wherein  the  mutation 
is  not  the  disease-initiating  event.  However,  due  to  the  pos¬ 
sibility  that  X  inactivation  in  some  of  our  female  patients 
may  be  skewed,  the  appearance  of  clonality  may  be  mis¬ 
leading.  Therefore,  this  data,  although  indicative  of 
JAK2V6,7F  not  being  the  initiating  event,  is  by  no  means 
sufficient  evidence  for  such  a  claim. 


Direct  evidence  for  JAK2V6nF  as  a  secondary  event 
Since  erythropoietin-independent  growth  of  BFU-E  colo¬ 
nies  (also  known  as  EEC)  is  considered  a  hallmark  of  PV 

[23] ,  we  decided  to  verify  if  our  assay  could  be  used  to  de¬ 
termine  the  JAK2  mutant  T-allele  frequency  in  single  colo¬ 
nies.  In  our  preliminary  experiments  using  JAK2 
sequencing  and  allele  discrimination  real-time  PCR,  we 
found  that  a  majority  of  untreated  PV  patients,  and  patients 
stable  on  their  respective  therapies,  had  erythropoietin-in¬ 
dependent  BFU-E  that  were  homozygous  for  the  JAK2V6,7F 
mutation.  Nevertheless,  some  colonies  were  heterozygous, 
while  rare  colonies  had  a  wild-type  JAK2  genotype  (data 
not  shown).  Using  our  novel  and  more  sensitive  quantitative 
assay  described  herein,  we  have  reexamined  this  issue. 
Since  individual  EEC  colonies  represent  the  “clonal” 
growth  and  proliferation  initiated  by  a  single  affected  cell 

[24] ,  allele  frequencies  should  reflect  those  normally  ob¬ 
served  for  any  germline  SNP  (0,  50,  and  100).  Mononuclear 
cells  from  two  female  and  two  male  PV  patients  were  used 
to  grow  and  harvest  89  individual  BFU-E  colonies  [25], 
a  total  of  69  genotypes  were  determined  (Table  3).  Devia¬ 
tion  from  the  expected  allele  frequency  was  small,  and  is 
hypothesized  to  be  due  to  contamination  by  other  colonies 
or  cells  picked  together  while  harvesting  the  colony  of  in¬ 
terest.  Colonies  with  allele  frequency  deviations  greater 
than  10%,  as  determined  by  the  JAK2  mutant  T-allele  fre¬ 
quency  distribution  (Supplemental  Fig.  S3),  were  discarded 
in  order  to  avoid  any  bias  during  genotyping  (Table  3).  Both 
heterozygous  and  homozygous  (mutant  and  wild-type) 
JAK2','6I7F  colonies  were  detected  (Table  3). 


Table  2.  Reproducibility  as  estimated  by  three  independent 
measurements 


Patient  number 

Date 

%T-allele  (PB-DNA)  %T-allele  (GNC-DNA) 

PV3 

4/13/2005 

16.4  (±  0.8) 

23.4 

PV4 

3/1/2005 

38.8  (±  3.0) 

60.4 

PV4 

6/30/2005 

16.8  (±  0.9) 

27.2  (±  4.7) 

PV4 

12/27/2005 

41.1  (±  2.8) 

55.1  (±  2.0) 

PV6 

8/31/2005 

N.D. 

0.0  (±  0.0) 

PV7 

1/4/2005 

74.8  (±  1.0) 

88.2  (±  0.7) 

PV7 

7/12/2005 

74.6  (±  0.6) 

84.8  (±  1.6) 

PV7 

10/14/2005 

67.3  (±  4.4) 

82.2  (±  1.6) 

PV8 

10/19/2005 

89.7  (±  0.3) 

97.4  (±  0.5) 

PV10 

5/1 1/2005 

63.4  (±  1.3) 

84.2  (±  1.0) 

PV12 

9/3/2005 

27.3  (±  1.9) 

40.7  (±  1.0) 

PV14 

1 1/30/2005 

74.7  (±  2.0) 

78.1  (±  1,1) 

PV15 

1  l/S/2005 

0.1  (±  0.0) 

0.0  (±  0.1) 

PV17 

12/8/2005 

37.1  (±  0.4) 

44.9  (±  1.3) 

PV18 

1 1/8/2005 

85.8  (±  0.8) 

91.4  (i  0.3) 

PV19 

11/15/2005 

66.5  (±  0.7) 

92.5  (±  4.9) 

PV21 

3/15/2005 

95.8  (±  0.4) 

99.9  (±  0.0) 

T-allele  frequency  (%)  was  determined  for  peripheral  blood  genomic  DNA 
(PB-DNA)  and  for  granulocyte  genomic  DNA  (GNC-DNA).  Results  are 
shown  as  the  mean  with  standard  deviation  in  between  parenthesis. 
N.D.,  not  determined. 
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%T  Expected  (LOG) 

Figure  2.  Limit  of  mutant  allele  detection  in  the  presence  of  a  large  excess 
of  wild-type  allele.  Genomic  DNA  from  HEL  cells  (JAK2  homozygous 
mutant)  was  mixed  in  decreasing  proportions  to  DNA  from  a  healthy  donor 
( JAK2  homozygous  wild-type).  Ensuing  mixtures  containing  10%,  1%, 
0.1%,  and  0.01%  HEL  gDNA  were  used  to  determine  JAK2V6,7F  mutant 
allele  frequency.  Linear  regression  of  kinetic  AS-PCR  determinations  vs 
expected  allele  frequencies  showed  excellent  correlation  (R2  =  0.996;  p 
<  0.002).  Data  represent  mean  ±  SD  of  three  independent  determinations. 


Discussion 

AS-PCR  allows  determination  of  the  allele  frequency  of 
any  polymorphism  or  mutation  based  on  the  difference  in 
the  number  of  PCR  cycles  needed  to  generate  detectable 
product.  AS-PCR  can  also  be  a  valuable  too!  to  determine 
the  frequency  of  rare  alleles  when  somatic  mutations  are 
present  in  a  small  percentage  of  cells  within  a  tissue  sam¬ 
ple.  Addition  of  mismatches  in  a  primer  requires  careful 
design,  strand  selection,  and  optimization  of  the  many  vari- 


Table  3,  JAK2V6I7F  allele  frequency  in  isolated  BFU-E  colonies 


Patient 

number 

T-allele  (%)  in  GNC 

BFU-E  with  0  U 
EPO 

BFU-E  with  3  U 
EPO 

G/G 

G/T 

T/T 

G/G 

G/T 

T/T 

PV18 

91.4 

0 

0 

7 

2 

1 

8 

PV22 

N.D. 

1 

1 

3 

5 

0 

7 

PV23 

20 

6 

1 

0 

7 

4 

0 

PV24 

89 

1 

2 

2 

6 

0 

4 

PV23 

55 

0 

0 

7 

0 

0 

5 

PV26 

42 

0 

0 

3 

2 

0 

1 

PV27 

31 

0 

2 

0 

2 

6 

0 

PV28 

91 

0 

0 

1 

0 

0 

0 

PV29 

61 

2 

0 

0 

0 

0 

0 

Total  colonies 

10 

6 

23 

24 

11 

25 

Genotypes  represent  the  quantitative  determination  of  the  frequency  of 
mutant  allele  with  cutoffs  as  follows:  GIG  (0-10%),  G/T  (40-60%),  and 
T/T  (90-100%).  Colonies  that  are  not  within  these  set  cutoffs  arc  deemed 
contaminated  by  other  cells  while  harvesting  and  are  discarded.  N.D.,  not 
determined. 


ables  that  affect  primer  stability  and  target  amplification. 
Furthermore,  the  use  of  LNA  in  the  design  of  oligonucleo¬ 
tide  primers  and  probes  resulted  in  improved  specificity  and 
allele  discrimination  compared  to  unmodified  primers 
[21,26].  However,  the  synergistic  effect,  as  shown  here,  of 
primers  harboring  a  mismatched  nucleoside  and  LNA  mod¬ 
ification  in  AS-PCR  has  not  been  investigated. 

We  have  demonstrated  here  that  use  of  AS-PCR  for 
quantitation  of  somatic  mutations  using  our  novel  primer 
design  is  sensitive,  accurate,  and  highly  reproducible.  Inde¬ 
pendent  investigators  have  tested  our  design  in  different 
laboratories  using  two  kinetic  PCR  platforms  (ABI-7000, 
made  by  Applied  Biosystems,  Foster  City,  CA,  USA,  and 
a  Light  Cycler,  made  by  Roche  Applied  Science,  Indianap¬ 
olis,  IN,  USA).  Results  for  the  samples  tested  were  nearly 
identical  between  platforms,  and  were  comparable  to  previ¬ 
ous  estimates  of  allele  frequency  obtained  by  pyrosequenc- 
ing  (results  not  shown).  Further,  amplicon  detection  during 
kinetic  AS-PCR  can  be  done  using  any  available  technol¬ 
ogy  (e.g.,  SYBR-Green,  TaqMan  Probe,  Molecular  Beacon, 
Hybridization  Probe)  at  the  discretion  and  convenience  of 
the  investigator. 

We  have  successfully  tested  our  primer  design  for  the 
quantitation  of  expressed  JAK2V6I7F  mutant  allele  mRNA 
from  PV  patient  samples  (data  not  shown).  Our  primer  de¬ 
sign  strategy  can  also  be  applied  to  determine  allelic  fre¬ 
quencies  of  other  SNP  or  mutations.  To  date,  we  have 
successfully  designed  five  different  primer/probe  sets  which 
will  be  used  to  detect  expression  of  X-chromosomal  SNPs 
that  are  subjected  to  inactivation  in  women,  for  determina¬ 
tion  of  clonality  [22] .  The  only  optimization  required,  when 
using  our  primer  design  strategy,  is  determination  of  the 
heterozygosity  correction  factor  for  proper  estimation  of  al¬ 
lele  frequencies. 

Potential  applications  for  our  primer  design  are  varied, 
such  as  quantitative  determination  of  mosaicism;  propor¬ 
tion  of  fetal  cells  in  maternal  circulation;  detection  of  min¬ 
imal  residual  disease  associated  with  known  somatic 
mutation  (such  as  reduction  of  malignant  cells  by  chemo¬ 
therapy  or  reappearance  of  resistant  clone);  rapid  monitor¬ 
ing  of  efficacy  of  new  drugs  in  “in  vitro”  systems  as  well  as 
clinical  trials;  and  many  others  that  require  quantitation  of 
allele  frequencies. 

We  used  our  quantitative  AS-PCR  assay  described  here 
to  investigate  the  order  of  genetic  events  leading  to  the 
JAK2V6’71'  mutation  in  the  ontogeny  of  clonal  hematopoie¬ 
sis  in  PV.  In  contrast  to  a  recent  study  reporting  the  absence 
of  73/C2vrt/7F-negative  EPO-independent  BFU-E  colonies 
[27],  we  found  such  colonies  in  samples  from  our  PV 
patients.  A  possible  explanation  for  this  discrepancy  may 
relate  to  the  sensitivity  and  methods  used  for  detection  of 
the  mutant  allele.  The  determination  of  zygosity  in  granulo¬ 
cytes  from  PV  patients  by  Scott  and  colleagues  was  done  in 
samples  that  were  collected  at  least  6  months  prior  to  the 
samples  used  for  clonogenic  assay  [27].  We  have  also  found 
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that  X-chromosome-based  methylation  clonality  assays 
have  significant  shortcomings  compared  to  transcrip¬ 
tional-based  X-chromosome  clonality  assays  [27],  In  addi¬ 
tion,  Scott  and  colleagues  used  a  qualitative  PCR  and 
restriction  endonuclease  digestion  assay  for  the  genotyping 
of  their  samples,  and  hence  could  not  distinguish  colonies 
potentially  contaminated  with  other  cells  [27].  In  this  assay, 
identification  of  JAK2  wild-type  genotypes  requires  the 
PCR  product  to  be  fully  digested  by  BsaXI.  Therefore,  ho¬ 
mozygous  JAK2  mutant  allele  and  heterozygous  colonies 
may  be  erroneously  genotyped  due  to  incomplete  digestion. 
Since  our  methods  make  use  of  real-time  monitoring  during 
PCR  amplification,  and  the  data  generated  is  quantitative, 
we  can  make  genotyping  determinations  with  greater  confi¬ 
dence.  Further,  the  in  vitro  analysis  of  BFU-E  response  to 
erythropoietin  (clonogenic  assay)  is  laborious  and  data 
not  always  easy  to  interpret;  however,  we  have  developed 
this  assay  and  have  routinely  used  it  since  1974  [23].  We 
have  validated  the  reproducibility  and  use  of  the  clonogenic 
assay  in  our  studies  of  a  large  number  of  congenital  and  ac¬ 
quired  polycythemic  states  such  as  primary  familial  and 
congenita]  polycythemia  and  Chuvash  polycythemia  [28]. 

The  presence  of  wild-type  JAK2  EPO-independent  colo¬ 
nies  suggests  the  existence  of  a  distinct,  and  as  of  now  un¬ 
known,  alternative  mechanism  active  in  erythropoiesis/ 
hematopoiesis.  Moreover,  our  results  indicate  the  presence 
of  an  undefined  molecular  lesion  that  precedes  the 
JAK2V6I7F  mutation.  Finally,  our  data  suggest  that  develop¬ 
ment  of  chemotherapeutics  that  target  the  JAK2V6!7F  clonal 
cells  may  not  suffice  to  cure  PV, 


Acknowledgments 

We  are  grateful  to  Yongli  Guan  for  technical  help  throughout  this 
work.  This  work  was  supported  by  a  1P01CA108671-O1A2 
MPD  Consortium  and  NIH  Grant  (HL-050077,  J.T.P.).  R.H.N. 
and  A.G.  were  supported  by  a  T-32  hematology  training  grant 
from  the  NIH  (DK-60445).  JJ.  was  supported  by  a  DOD  Grant 
W81XWH-05- 1-0535. 


References 

1.  Skoda  R,  Prchal  JT.  Chronic  myeloproliferative  disorders — introduc¬ 
tion,  Semin  Hematol.  2005;42:181-183. 

2.  Baxter  EJ,  Scott  LM.  Campbell  PJ,  et  al.  Acquired  mutation  of  the  ty¬ 
rosine  kinase  JAK2  in  human  myeloproliferative  disorders.  Lancet. 
2005;365:1054-1061. 

3.  James  C,  Ugo  V,  Lc  Couedic  JP,  et  al.  A  unique  clonal  JAK2  mutation 
leading  to  constitutive  signalling  causes  polycythaemia  vera.  Nature. 
2005;434:1144-1 148. 

4.  Krnlovics  R,  Passamonti  F,  Buser  AS,  et  al.  A  gain-of-function  muta¬ 
tion  of  JAK2  in  myeloproliferative  disorders.  N  Engl  J  Med.  2005;352: 
1779-1790. 

5.  Levine  RL,  Wadleigh  M,  Cools  J,  et  al.  Activating  mutation  in  the  ty¬ 
rosine  kinase  JAK2  in  polycythemia  vera,  essential  thrombocythemia, 
and  myeloid  metaplasia  with  myelofibrosis.  Cancer  Cell.  2005;7:387- 
397. 


6.  Zhao  R,  Xing  S,  Li  Z,  et  al.  Identification  of  an  acquired  JAK2 
mutation  in  polycythemia  vera.  J  Bioi  Chem.  2005;280:22788-22792. 

7.  Ugo  V,  Marzac  C,  Teyssandier  I,  et  al.  Multiple  signaling  pathways  are 
involved  in  erythropoietin-independent  differentiation  of  erythroid 
progenitors  in  polycythemia  vera.  Exp  Hematol.  2004;32:179-187. 

8.  Lacout  C,  Pisani  DF,  Tulliez  M,  Moreau  Gachelin  F,  Vainchcnkcr  W, 
Villeval  JL.  JAK2V617F  expression  in  murine  hematopoietic  cells 
leads  to  MPD  mimicking  human  PV  with  secondary  myelofibrosis. 
Blood.  2006;108:1652-1660. 

9.  Kralovics  R,  Guan  Y,  Prchal  JT.  Acquired  uniparental  disomy  of  chro¬ 
mosome  9p  is  a  frequent  stem  cell  defect  in  polycythemia  vera,  Exp 
Hematol.  2002;30:229-236. 

10.  Kralovics  R,  Stockton  DW,  Prchal  JT.  Clonal  hematopoiesis  in  famil¬ 
ial  polycythemia  vera  suggests  the  involvement  of  multiple  mutational 
events  in  the  early  pathogenesis  of  the  disease.  Blood.  2003;  102:3793- 
3796, 

1 1 .  Kralovics  R,  Teo  SS,  Li  S,  et  al.  Acquisition  of  the  V61 7F  mutation  of 
JAK2  is  a  late  genetic  event  in  a  subset  of  patients  with  myeloprolif¬ 
erative  disorders.  Blood.  2006;108:1377-1380. 

12.  Bellanne-Chantelot  C,  Chaumarel  I,  Labopin  M,  et  al.  Genetic  and 
clinical  implications  of  the  Val617Phe  JAK2  mutation  in  72  families 
with  myeloproliferative  disorders.  Blood.  2006;108:346-352. 

13.  Germer  S,  Holland  MJ,  Higuchi  R.  High-throughput  SNP  allele-fre¬ 
quency  determination  in  pooled  DNA  samples  by  kinetic  PCR.  Ge¬ 
nome  Res.  2000;10:258-266. 

14.  Jelinek  J,  Prchal  JT.  Oxygen-dependent  regulation  of  erythropoiesis. 
Methods  Enzymol.  2004;381:201-210. 

15.  Newton  CR,  Graham  A,  Heptinstall  LE,  et  al.  Analysis  of  any  point 
mutation  in  DNA.  The  amplification  refractory  mutation  system 
(ARMS).  Nucleic  Acids  Res.  1989;17:2503-2516. 

16.  Petersen  M,  Nielsen  CB,  Nielsen  KE,  et  al.  The  conformations  of 
locked  nucleic  acids  (LNA).  J  Mol  Recognit.  2000;13:44-53. 

17.  Jelinek  J1,  Oki  Y,  Gbaribyan  V,  et  al.  JAK2  mutation  1849G>T  is  rare 
in  acute  leukemias  but  can  be  found  in  CMML,  Philadelphia  chromo- 
some-negative  CML,  and  megakaryocylic  leukemia.  Blood.  2005; 
106:3370-3373. 

18.  Ahtnadian  A,  Gharizadeh  B,  O’Meara  D,  Odebcrg  J,  Lundeberg  J. 
Genotyping  by  apyrase-mediated  allele-specific  extension.  Nucleic 
Acids  Res.  2001;29:E!21. 

19.  Bottema  CD,  Sommer  SS.  PCR  amplification  of  specific  alleles:  rapid 
detection  of  known  mutations  and  polymorphisms,  Mutat  Res.  1993; 
288:93-102. 

20.  Kaltenbock  B,  Schneider  R.  Differential  amplification  kinetics  for  point 
mutation  analysis  by  PCR.  Biotechniques.  1998;24:202-204, 206. 

21.  Latorra  D,  Campbell  K,  Wo  Her  A,  Hurley  JM.  Enhanced  allele-spe¬ 
cific  PCR  discrimination  in  SNP  genotyping  using  3'  locked  nucleic 
acid  (LNA)  primers.  Hum  Mutat.  2003;22:79-85. 

22.  Liu  E,  Jelinek  J,  Pastore  YD,  Guan  Y,  Prchal  JF,  Prchal  JT.  Discrim¬ 
ination  of  polycythemias  and  thrombocytoses  by  novel,  simple,  accu¬ 
rate  clonality  assays  and  comparison  with  PRV-1  expression  and 
BFU-E  response  to  erythropoietin.  Blood.  2003;101:3294-3301. 

23.  Prchal  JF,  Axelrad  AA.  Letter:  Bone-marrow  responses  in  polycythe¬ 
mia  vera.  N  Engl  J  Med.  1974;290: 1 382. 

24.  Prchal  JF,  Adamson  JW,  Steinmann  L,  Fialkow  PJ.  Human  erythroid 
colony  formation  in  vitro:  evidence  for  clonai  origin.  J  Cell  Physiol. 
1976;89:489^t92, 

25.  Stopka  T,  Zivny  JH,  Stopkova  P,  Prchal  JF,  Prchal  JT.  Human  hemato¬ 
poietic  progenitors  express  etythropoietin.  Blood.  1998;91:3766-3772. 

26.  Latorra  D,  Arar  K,  Hurley  JM.  Design  considerations  and  effects  of 
LNA  in  PCR  primers.  Mol  Cell  Probes.  2003:17:253-259. 

27.  Scott  LM,  Scott  MA,  Campbell  PJ,  Green  AR.  Progenitors  homozy¬ 
gous  for  the  V617F  mutation  occur  in  most  patients  with  polycythemia 
vera,  but  not  essential  thrombocythemia.  Blood.  2006;108:2435-2437. 

28.  Prchal  JT.  Polycythemia  vera  and  other  primary  polycythemias.  Curr 
Opin  Hematol.  2005;12:112-116. 


Page  35 


R.H.  Nussenzveig  et  aU  Experimental  Hematology  35  (2007)  32-38 


38.el 


Table  SI.  Primer  and  Probe  sequences 

(A)  Primer/Probe  Modification  <B)  Sequence  5'  to  3' 


FAM-AS-JAK2-MGB 
gJAK2-F 
G-gJAK2-R 
T-  gJAK2-R 

G-gJAK2-LNA-(no- 1  mut)-R 
T-  gJAK2-LNA-(no-lmut)-R 
G-gJAK2-(— lmut)-R 
T-  gJAK2-(— lmut)-R 
G-gJAK2-{0)LNA-R 
T-  gJAK2-(0)LNA-R 
G-gJAK2-(—  1  )LNA-R 
T-  gJAK2-(—  1  )LNA-R 
G-gJAK2-(— 2)LNA-R 
T-  gJAK2-(-2)LNA-R 
G-gJAK2-(-3)LNA-R 
T-  gJAK2-(— 3)LNA-R 
G-gJAK2-(-4)LNA-R 
T-  gJAK2-(— 4)LNA-R 


FAM  labeled  MGB  probe  with  non-fiuorcscent  quencher  (NFQ) 

N/A 

None 

None 

LNA  at  -2 
LNA  at  -2 

—  1  Mismatch 

—  1  Mismatch 

LNA  at  0  and  —  1  Mismatch 
LNA  at  0  and  -1  Mismatch 
LNA  at  -l  and  —1  Mismatch 
LNA  at  -1  and  -1  Mismatch 
LNA  at  -2  and  —1  Mismatch 
LNA  at  —2  and  —1  Mismatch 
LNA  at  -3  and  —1  Mismatch 
LNA  at  —3  and  - 1  Mismatch 
LNA  at  —4  and  —1  Mismatch 
LNA  at  -4  and  -1  Mismatch 


5'  6FAM  -  CTTGCTCATCATACTTGC  -  MGBNFQ  3' 
TTATGGACAACAGTCAAACAACAAT 
TTTACTTACTCTCGTCTCCACAGAC 
TTTACTTACTCTCGTCTCCACAGAA 
TTTACTTACTCTCGTCTCCACA  G  AC 
TTTACTTACTCTCGTCTCCACA  G  AA 
TTTACTTACTCTCGTCTCCACAG  t  C 
TTTACTTACTCTCGTCTCCACAG  t  A 
TTTACTTACTCTCGTCTCCACAG  t  C 
TTTACTTACTCTCGTCTCCACAG  t  A 
TTTACTTACTCTCGTCTCCACAG  t  C 
TTTACTTACTCTCGTCTCCACAG  /  A 
TTTACTTACTCTCGTCTCCACA  (ft C 
TTTACTTACTCTCGTCTCCACA  G  t  A 
TITACTTACTCTCGTCTCCAC  A  G  t  C 
TTTACTTACTCTCGTCTCCAC  AGtA 
TTTACTTACTCTCGTCTCCA  C  AG  t  C 
TTTACTTACTCTCGTCTCCA  C  AGtA 


A)  Allele  specific  primers  are  named  beginning  with  uppercase  G  or  T  based  on  the  sequence  of  the  sense  strand.  The  G-primers  are  specific  for  the  wild-type 
allele,  while,  the  T-primers  are  specific  for  the  mutant  allele. 

B)  LNA  bases  arc  depicted  in  uppercase  italics,  bolded  and  underlined;  while,  the  —1  mismatch  is  highlighted  in  bold  lowercase.  Base  positions  are  counted 
beginning  at  the  3’  end  of  the  primers. 


5‘ - t-W-C  3' 

OR 

S' - L-H-A  3' 


G 


Figure  SI.  Diagram  of  the  mechanism  of  action  for  our  primer  design.  The  diagram  depicts  the  strand  used  for  allele-specific  priming  and  extension.  The 
strand  primed  and  extended  by  the  universal  primer  is  not  shown.  In  the  first  step,  two  different  reactions  arc  prepared,  each  containing  one  of  the  allele- 
specific  primers  and  the  universal  primer.  In  the  second  step,  after  heat  denaturation,  the  allele-specific  primer  anneals  and,  depending  on  3’  end  hybridiza¬ 
tion,  Taq  DNA  polymerase  extension  occurs.  Finally,  in  the  third  step,  once  extension  takes  place,  the  5’  to  3’  exonuclease  activity  of  Taq  cleaves  the  probe 
between  the  reporter  and  quencher.  Subsequent  PCR  cycles  favor  the  newly  synthesized  strand  since  it  incorporates  the  -1  mismatched  nucleoside.  We  hy¬ 
pothesize  that  the  LNA,  in  our  primer  design,  helps  stabilize  3’  end  annealing  by  improving  base  pair  stacking  of  matched  bases,  increasing  primer  7m,  and 
decreasing  AG.  Hence,  as  depicted  above,  3'  end  floppiness,  resulting  from  the  introduction  of  the  —1  mismatch,  is  compensated  by  the  presence  of  the  LNA, 
leading  to  enhanced  reproducibility,  specificity  and  efficiency  of  the  reaction. 
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Pyrosequencing  (%T -allele) 


Figure  S2.  Linear  regression  of  the  correlation  between  AS-PCR  vs. 
Pyrosequencing.  JAK2V6I7F  (%T-allele)  frequency,  for  the  same  patient 
sample,  were  independently  determined  by  real-time  AS-PCR  and  pyrose¬ 
quencing.  An  excellent  correlation  between  both  assay  platforms  was 
found  in  a  linear  correlation  model  as  shown  above. 


Percent  of  T-allele  in  individual  colonies 


Figure  S3.  BFU-E  colony  analysis  showing  a  trimodal  distribution  of 
T-allele  frequencies.  Frequency  distribution  of  all  BFU-E  colonies  picked 
and  used  for  genotyping  the  JAK2V6/7F  mutation,  as  described  in  the 
accompanying  manuscript.  The  cutoff's  used  for  defining  genotypes  are 
indicated  in  table  3. 
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Objective.  Polycythemia  rubra  vera-1  (PRV-1)  is  a  GPI-Iinked  protein  that  is  expressed  on 
a  subgroup  of  neutrophils.  The  number  of  PRV- 1-expressing  neutrophils  increases  in  preg¬ 
nancy  and  sepsis,  or  after  administration  of  granulocyte  colony-stimulating  factor.  Expression 
of  the  PRV-1  gene  is  also  increased  in  patients  with  polycythemia  vcra  (PV)  and  essential 
thrombocythemia  (ET).  We  investigated  whether  DNA  methylation  of  the  PRV-1  gene  has 
a  role  in  regulation  of  transcription  and  expression  of  the  PRV-1  protein. 

Methods.  We  compared  the  level  of  methylation  of  the  PRV-1  gene  and  expression  of  the  PRV- 
1  mRNA  in  normal  neutrophils  expressing  PRV-1  to  those  that  are  PRV-l-negative.  We  also 
studied  PRV-1  methylation  and  mRNA  expression  in  patients  with  Philadelphia  chromosome- 
negative  myeloproliferative  disorders  and  in  an  in  vitro  model  of  DNA  demethylation. 
Results.  We  found  that  methylation  of  CpG  dinucleotides  close  to  initiation  codon  of  the  PRV- 
1  gene  was  inversely  related  to  expression  of  PRV-1  in  normal  neutrophils.  Furthermore,  over¬ 
expression  of  the  PRV-1  gene  in  PV  and  ET  is  associated  with  a  decrease  in  methylation  of  this 
gene.  Among  patients  with  PV  and  ET,  methylation  of  the  PRV-1  gene  is  also  inversely  corre¬ 
lated  with  the  presence  of  the  JAK2V617F  somatic  mutation.  In  an  in  vitro  model,  exposure  of 
KG1  and  KGla  cells  to  a  DNA  demcthylating  agent  caused  a  decrease  in  methylation  of  the 
PRV-1  gene  and  increased  its  mRNA  level. 

Conclusion.  DNA  methylation  regulates  PRV-1  expression  under  physiologic  and  pathologic 
conditions.  ©  2007  1SEH  -  Society  for  Hematology  and  Stem  Cells.  Published  by  Elsevier 
Inc. 


Polycythemia  rubra  vera-1  (PRV-1),  NB1,  or  CD  177  is 
a  GPI-linked  surface  protein  that  is  expressed  on  neutro¬ 
phils  [1J.  Interestingly,  PRV-1  is  normally  expressed  only 
on  a  subgroup  of  neutrophils  [2,3].  Circulating  neutrophils 
can  be  divided  into  PRV- 1-positive  and  PRV-l-negative 
subgroups  according  to  the  presence  or  absence  of  PRV-1 
on  their  surface  (Fig.  1 ).  Regulation  of  this  unique  expres¬ 
sion  pattern  is  at  the  level  of  transcription  of  the  PRV- 1  gene, 
evident  by  absence  of  PRV-1  mRNA  in  PRV-l-negative 
neutrophils  [3].  The  relative  proportion  of  PRV-1 -positive 
neutrophils  varies  from  one  individual  to  another,  and 
even  in  the  same  individual  over  time  and  under  different 
physiologic  and  pathologic  stimuli.  Sepsis,  pregnancy, 
and  administration  of  granulocyte  colony-stimulating  factor 
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increase  expression  of  PRV- 1  [3-6].  More  recently,  it  was 
recognized  that  PRV-1  mRNA  level  is  increased  substan¬ 
tially  in  Philadelphia  chromosome-negative  myeloprolifer¬ 
ative  disorders,  especially  in  polycythemia  vera  (PV)  [7]. 
We  think  that  mechanisms  regulating  PRV-1  expression 
might  also  be  involved  in  pathogenesis  of  myeloprolifera¬ 
tive  disorders;  and  studying  these  mechanisms  might  be 
helpful  for  a  better  understanding  of  the  molecular  changes 
of  myeloproliferative  disorders.  Additionally,  we  and  others 
have  shown  that  PRV-1  regulates  cell  proliferation  [5,8,9], 
and  might  play  a  direct  role  in  pathogenesis  of  myeloprolif¬ 
erative  disorders.  We  investigated  whether  methylation  of 
the  PRV-1  gene  affects  its  transcription  and  has  any  role 
in  biphenotypic  expression  pattern  of  PRV-1  in  neutrophils. 
The  PRV- 1  gene  is  located  on  chromosome  19  and  is  devoid 
of  CpG  islands.  We  studied  cytosine  methylation  of  CpG 
dinucleotides  close  to  the  PRV-1  gene  transcription  start 
site  at  cytosines  30,  41,  and  43  (six  nucleotides  before, 
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FITC  (PRV-1 ) 

Figure  1.  Polycythemia  rubra  vern-1  (PRV-1)  expression  on  normal  neu¬ 
trophil.  Results  of  flow  cytometry  on  neutrophils  from  a  normal  donor,  us¬ 
ing  MEM  166  (anti  PRV-1  antibody)  shown  as  a  frequency  distribution 
plot.  Neutrophils  had  a  biphenotypic  expression  pattern  for  PRV-1,  and 
in  this  donor  71%  of  all  neutrophils  expressed  PRV-1.  FITC  =  fluorescein 
isothiocyunate. 

three  and  five  nucleotides  after,  the  initiation  codon,  respec¬ 
tively)  (Fig.  2).  We  compared  CpG  methylation  between 
PRV-1 -positive  and  negative  neutrophils  in  a  group  of 
healthy  subjects.  We  further  studied  PV  as  a  mode!  for  dys- 
regulation  of  the  PRV-1  gene  expression.  We  investigated 
the  correlation  between  the  PRV-1  mRNA  level  and  CpG 
methylation  of  its  gene  in  neutrophils  of  patients  with  PV, 
essential  thrombocythemia  (ET),  primary  familial  or  con¬ 
genital  polycythemia;  and  in  a  group  of  normal  individuals. 
We  found  that  a  lower  methylation  of  CpG  sites  close  to 
transcription  start  site  was  associated  with  higher  PRV-1 
expression  both  in  normal  individuals  and  in  patients  with 
PV  and  ET. 


Materials  and  methods 

Isolation  of  neutrophil  subgroups 
and  determination  of  the  PRV-I  mRNA  levels 
To  compare  DNA  methylation  among  different  neutrophil  PRV-1 
subgroups  in  normal  individuals,  blood  samples  from  10  healthy 
subjects  were  collected  into  3.5%  citrate  anticoagulant  at  a  1:6 
ratio.  Using  Histopaque  (Sigma,  St  Louis,  MO,  USA)  density- 
gradient  centrifugation  followed  by  ammonium  chloride/bicarbon¬ 
ate  lysis  of  red  blood  cells,  granulocytes  were  separated  from  the 
other  blood  cells  [10],  To  separate  neutrophil  granulocytes  into  dif¬ 
ferent  PRV-  1-expression  subgroups,  they  were  labeled  with  fluores¬ 
cein  isothiocyanate-conjugated  anti-PRV-1  antibody  (MEM166, 
Serotec)  and  phycoerythrin-conjugated  anti-CD  lib  antibody 
(Dako,  Carpinteria,  CA,  USA).  Using  cell-sorter  EPICS  ALTRA 
(Beckman  Coulter,  Fullerton,  CA,  USA),  CDllb/CD18  (MAC- 
l)-expressing  cells  in  the  size  range  of  granulocytes  were  separated 
into  PRV-l-positive  and  PRV- 1 -negative  subgroups  (Fig.  1),  We 
used  Qiagen  blood  DNA  kit  to  extract  DNA  from  neutrophils 
(106  cells/sample)  and  used  bisulfite  treatment  to  prepare  the 
DNA  samples  for  methylation  studies  [1 1], 

To  study  the  correlation  between  the  extent  of  PRV-1  methyl¬ 
ation  and  its  expression  in  patients  with  PV,  we  used  DNA  and 
RNA  samples  from  unsorted  granulocytes  of  19  patients  with 
PV,  12  with  ET,  and  1 1  with  primary  familial  and  congenital  poly¬ 
cythemia.  DNA  samples  were  prepared  for  methylation  studies 
and  RNA  samples  were  used  for  determining  the  PRV-1  mRNA, 
as  described  previously  [11,12],  All  of  the  blood  samples  were  ob¬ 
tained  according  to  protocol  approved  by  the  institutional  review 
board  and  after  subjects  signed  an  informed  consent. 

In  vitro  exposure  of  cells  to  DNA  demethylating  agent 
Leukemia  cell  lines  KG-1  and  KG-Ia  were  obtained  from  the 
American  Type  Culture  Collection  (Manassas,  VA,  USA).  Cells 
were  grown  in  Iscove’s  modified  Dulbecco's  medium  supple¬ 
mented  with  fetal  bovine  serum  (20%  v/v).  DNA  demethylation 
was  induced  by  a  4-day  treatment  with  200  nM  2-deoxy-5-azacy- 
tidine  [I  I]  (DAC;  Sigma).  Fresh  DAC  was  added  to  exponentially 


j  CpG  dinucleotide 
r=i  Exon 


Figure  2.  CpG  dinucleotides  position  in  relation  to  transcription  start  site  of  the  polycythemia  rubra  vera-I  (PRV-1 )  gene.  Exons  are  indicated  as  boxes  and 
CpG  sites  as  vertical  bars.  Three  CpG  dinucleotides  are  located  around  the  initiation  codon  in  the  exon  1  of  the  PRV-1  gene:  cytosine  (C)‘10pG,  C^'pG,  and 
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growing  cells  at  days  0,  1,  2,  and  3.  Parallel  cultures  of  mock- 
treated  cells  were  kept  as  controls.  Cells  were  harvested  at  day 
4,  lysed  in  Trizol  (Invitrogen,  Carlsbad,  CA,  USA),  and  RNA 
plus  DNA  were  isolated. 

DNA  sequencing  and  methylation  studies 

We  determined  the  DNA  sequence  of  a  fragment  of  genomic  DNA 
encompassing  the  initiation  codon  of  the  PRV-1  gene  to  determine 
the  genotype  at  the  known  polymorphism  at  guanosine  42  (G42C) 
[13]  and  thus  the  presence  or  absence  of  the  second  CpG  site  at 
C41.  After  bisulfite  treatment,  the  same  region  of  the  genomic 
DNA  was  amplified  and  cloned  into  pCR4  sequencing  vector  (In¬ 
vitrogen).  For  each  polymerase  chain  reaction  (PCR)  product,  25 
clones  representing  individual  alleles  were  pyrosequenced  to 
determine  the  methylation  status  of  cytosine  residues  at  three 
CpG  sites;  CpG  30,  41,  and  43,  close  to  the  PRV-1  transcription 
start  site  (numbered  according  to  the  sequence  accession  no. 
AF146747  from  the  GenBank  data)  (Fig.  2)  [11].  To  examine 
methylation  status  of  CpG  30  in  patient  samples,  bisulfite  PCR 
products  were  directly  analyzed  by  pyrosequencing.  Sequences 
of  all  primers  are  available  in  the  supplementary  data  (Supplemen¬ 
tary  Table).  We  also  determined  mutation  status  of  the  Janus 
kinase  2  ( JAK2 )  gene  at  codon  617  and  quantified  the  mutant 
allele  by  pyrosequencing  as  described  previously  [14]. 

Statistical  analysis 

We  used  Fisher’s  exact  test,  Chi-square  test,  Mann-Whitney  non- 
parametric  test,  and  nonparamelric  Spearman  correlation  for  data 
analysis. 


Results 

Expression  of  PRV-1  on  normal 
neutrophils  and  methylation  of  the  PRV-1  gene 
We  compared  methylation  of  the  PRV-1  gene  in  pairs  of 
PRV- 1-positive  and  PRV- 1-negative  neutrophil  subpopula¬ 
tions  isolated  from  10  normal  individuals.  Neutrophils 
expressing  PRV-1  surface  antigen  showed  consistently 
lower  methylation  levels  at  all  CpG  sites  near  the  PRV-1 
transcription  start  site  compared  to  PRV-negative  ones 
(Fig.  3).  We  sequenced  509  PCR  clones  (25  clones  per  sam¬ 
ple)  to  analyze  methylation  status  of  CpG  dinucleotides  on 
individual  alleles.  On  average,  C30  of  the  PRV-1  gene  was 
13%  ±  2%  methylated  in  PRV-1 -positive  neutrophils  and 
33%  ±  4%  in  PRV- 1-negative  neutrophils  (p  —  0.0015, 
Mann-Whitney  test).  Average  methylation  levels  of  C41  ad¬ 
jacent  to  C42G  polymorphic  site  were  52%  ±  7%  in  PRV- 
1-positive  and  71%  ±  3%  in  PRV- 1-negative  neutrophils 
(p  —  0.0041).  C43  showed  average  methylation  21%  ±  4% 
in  PRV-1 -positive  and  33%  ±  3%  in  PRV-1 -negative  neu¬ 
trophils;  this  difference  was  not  significant.  Cloning  and  se¬ 
quencing  approach  enabled  to  detect  methylation  pattern  of 
individual  alleles.  Frequencies  of  alleles  showing  all  possi¬ 
ble  combinations  of  methylation  status  of  three  examined 
CpG  sites  are  shown  in  Figure  4.  There  was  a  statistically 
significant  predominance  of  completely  unmethylated 


Figure  3.  Comparison  of  methylation  of  CpG  dinucleotides  between  poly¬ 
cythemia  rubra  vera-1  (PRV- 1  )-positive  and  -negative  neutrophils.  The 
methylation  status  of  CpG  dinucleotides  at  cytosine  (C)30,  C41,  and 
C43  sites  downstream  from  the  transcription  start  site  of  the  PRV-1  gene 
was  compared  between  PRV-1 -positive  (full  circles)  and -negative  neutro¬ 
phils  (open  circles)  in  each  subject  (total  of  10  normal  control  subjects). 
Asterisks  denote  significant  differences  (Mann-Whitney  nonparametric 
test,  p  <  0.05). 

alleles  in  PRV- 1-positive  granulocytes  while  alleles  with 
the  first  two  sites  methylated  or  with  all  three  sites  methyl¬ 
ated  were  predominant  in  PRV-1 -negative  granulocytes. 

CpG  methylation  and  the  G42C 
polymorphism  of  the  PRV-1  gene 

One  of  the  single  nucleotide  polymorphisms  in  the  PRV-1 
gene  is  due  to  replacement  of  guanosine  42  (major  allele 
with  frequency  of  54-68%)  by  a  cytosine  (minor  allele 
with  a  frequency  of  32-46%)  [6,13],  We  analyzed  methyl¬ 
ation  pattern  of  individual  42G  or  42C  alleles  based  on  the 
genotype  at  the  G42C  polymorphic  site.  There  was  no  sig¬ 
nificant  difference  in  C30  methylation  between  the  geno¬ 
types:  25  of  150  methylated  alleles  (17%)  for  the  42C/C 
genotype,  51  of  174  methylated  alleles  (23%)  for  42C/G, 
and  37  of  134  methylated  alleles  (28%)  for  the  42G/G  ge¬ 
notype.  Striking  differences  depending  on  G42C  genotype 
were  seen  in  C43  methylation.  Methylation  at  C43  was  ob¬ 
served  in  8  of  150  (8%)  alleles  in  granulocytes  with  the 
42C/C  genotype,  in  29  of  225  (29%)  alleles  in  granulocytes 
with  the  42C/G  genotype,  and  in  42  of  134  (42%)  alleles  in 
granulocytes  with  the  42G/G  genotype  (p  <  0.0001 ,  Chi- 
square  test).  To  avoid  the  effect  of  the  G42C  polymorphism 
on  our  methylation  analysis,  we  focused  on  the  methylation 
status  of  C30  site  in  subsequent  studies. 

Methylation  of  PRV-1  in  patients  with  PV,  ET, 
or  primary  familial  and  congenital  polycythemia 
Having  established  the  importance  of  C30  methylation  in 
sorted  neutrophils  from  normal  controls,  we  examined 
C30  methylation  and  PRV-1  mRNA  levels  in  samples  of 
unsorted  neutrophils  from  patients  with  PV,  ET,  primary  fa¬ 
milial  or  congenital  polycythemia,  and  in  normal  subjects. 
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Figure  4.  Methylation  pattern  of  CpG  sites  cytosine  (C)  30,  C41,  and  C43 
in  individual  alleles  of  the  polycythemia  rubra  vera-1  (PRV-1)  gene.  The 
methylation  status  of  the  PRV-1  gene  was  determined  in  PRV- 1-positive 
and  -negative  granulocytes  of  10  normal  subjects.  Black  bars  denote  allelic 
frequency  in  PRV-1 -positive  granulocytes,  white  bars  in  PRV-1 -negative 
granulocytes.  U  =  unmethylated,  M  methylated  CpG  site.  Missing  C41 
CpG  sites  in  42C/C  homozygotes  arc  indicated  by  an  underscore  sign. 
^Statistically  significant  differences  (p  <  0.05,  Fisher’s  exact  test). 

We  found  that  the  PRV-1  mRNA  level  was  inversely  corre¬ 
lated  to  the  extent  of  methylation  of  the  PRV-1  gene  at  C30 
in  PV  (Spearman  r  —  —0.86,  p  <  0.0001)  and  ET  patients 
(r  =  —0.64,  p  —  0.02)  (Fig.  5).  Patients  with  overexpres¬ 
sion  of  PRV-1  had  a  lower  extent  of  methylation  of  the 


PRV-1  gene  at  C30  ( <  20%  methylated  cytosine)  compared 
to  patients  with  normal  expression  of  PRV-1  (40-60% 
methylation). 

Methylation  of  PRV-1  gene  and  the  JAK2V6,7F  mutation 
Because  PRV-1  mRNA  is  overexpressed  in  PV  and  ET  pa¬ 
tients  carrying  JAK2V617F  mutation,  we  also  examined  the 
relationship  between  JAK2V6,7F  mutational  status  and  C30 
methylation  of  the  PRV-l  gene.  Samples  from  patients  with 
JAK2V6I7F  mutation  showed  significantly  lower  methylation 
of  C30  site,  when  compared  to  individuals  with  normal  JAK2 
sequence  (Fig.  6),  suggesting  that  constitutive  activation  of 
JAK2  signaling  results  in  PRV-1  demethylation.  Methylation 
of  PRV-1  at  C30  inversely  correlated  with  PRV-1  expression 
in  16  patients  harboring  5%  to  92%  JAK2617F  mutant  allele 
(Spearman  r  —  -0.78,  p  =  0.0004).  The  burden  of 
JAK26J7F  mutant  allele  was  also  inversely  related  to  C30 
methylation  of  PRV-1,  however,  the  correlation  was  not 
statistically  significant  (Spearman  r  —  —0.28,  p  ~  0.3). 

In  vitro  model  for  demethylation  of  the  PRV-1 
gene  and  its  effect  on  expression  of  PRV-1  mRNA 
We  studied  DNA  methylation  and  mRNA  expression  of 
PRV-1  in  human  leukemia  cell  lines  (KG1  and  KGla)  be¬ 
fore  and  after  4  days  of  exposure  to  a  demethylating  agent, 
2-deoxy-5-azacytidine  (decitabine  or  DAC),  in  the  cell  cul¬ 
ture  media.  ICG-1  and  ICG-la  myeloid  leukemia  cell  lines 
showed  high  methylation  levels  of  C30  and  low  levels  of 
PRV-1  mRNA  expression.  Treatment  with  decitabine,  de¬ 
creased  methylation  level  of  the  PRV-1  gene  in  KG1  and 
KGla  cells  (35%  and  44%,  respectively)  and  this  was  asso¬ 
ciated  with  a  drastic  increase  in  the  amount  of  PRV-1 
mRNA  in  these  cells  (28-fold  and  17-fold  increases  in 
mRNA  concentration,  respectively)  (Table  1). 


Figure  5.  Methylation  of  the  polycythemia  rubra  vera-1  (PRV-1)  gene  and  expression  of  PRV-1  in  patients  with  polycythemia  vera  (PV)  and  essential  throm- 
bocylhemia  (ET).  The  relative  amount  of  PRV-1  mRNA  to  a  control  mRNA  in  neutrophils  of  patients  with  PV  (n  =  19)  or  ET  (n  =  12)  was  determined  by 
real-time  polymerase  chain  reaction  and  is  shown  in  the  Y-axis.  The  percentage  of  methylation  of  C30  in  neutrophils  of  the  same  patients  was  determined  by 
pyrosequencing  and  is  shown  in  the  X-axis.  Solid  line  represents  linear  regression  and  broken  lines  show  95%  confidence  intervals.  Negative  correlation 
between  PRV-1  mRNA  levels  and  C30  methylation  was  observed  in  PV  patients,  p  <  0.0001  (Spearman  nonparametric  correlation  coefficient  was 
-0.86)  and  in  ET  patients,  p  =  0.02  (Spearman  nonparametric  correlation  coefficient  was  -0.64). 


J.  Jelinek  et  al.l  Experimental  Hematology  35  (2007)  1677— 1683 


1681 


Figure  6.  Methylation  of  the  polycythemia  rubra  vera-1  ( PRV-1 )  gene  and 
the  JAK2v6I7F  mutation.  The  percentage  of  methylation  of  cytosine  (C)  30 
in  the  PRV-1  gene  was  determined  in  subjects  with  and  without  the 
JAK2w''n  mutation.  The  group  homozygous  for  wild-type  JAK2  allele  in¬ 
cluded  40  blond  samples  (15  controls,  11  primary  familial  and  congenital 
polycythemia,  9  essential  thrombocythemia  [ET],  and  4  polycythemia  vera 
[PV])  and  the  JAK2V6'7F  group  included  18  blood  samples  (3  ET  and  15 
PV).  Nonparametric  Mann-Whitney  test  was  used  to  compare  the  two 
groups  (p  =  0.02). 

Discussion 

Presence  of  a  biphenotypic  expression  pattern  of  PRV-1  on 
circulating  neutrophils  cannot  be  readily  explained  by  the 
presence  of  DNA  polymorphisms,  different  stage  of  matu¬ 
ration  of  neutrophils,  or  transcription  factor  profile.  Circu¬ 
lating  neutrophils  are  cells  derived  from  the  same  tissue 
of  origin,  have  the  identical  stage  of  maturation,  and  are 
controlled  by  identical  transcription  factors  and  chemo- 
kines.  We  hypothesized  that  DNA  methylation  may  have 
a  role  in  bimodal  expression  of  the  PRV-1  gene  in  neutro¬ 
phils,  DNA  methylation  is  an  important  epigenetic  mecha¬ 
nism  regulating  gene  expression  during  development  and 
maturation  of  an  organism.  Usually,  methylation  status  of 
cytosine  base  of  CpG  dinucleotides  in  CpG-rich  islands 
of  a  gene  promoter  affects  transcription  of  that  gene.  How¬ 
ever,  there  are  several  examples  of  cytokine  gene  regulation 


Tbbic  1.  Decitabine  treatment  caused  demethylation  and  induced 
expression  of  PRV- 1  gene 
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Polycythemia  rubra  vera  —  1  (PRV-1)  mRNA  expression  was  measured  as 
ACt  -  the  difference  of  cycles  at  threshold  0.2  for  glyccraldehyde  phos¬ 
phate  dehydrogenase  and  PRV- 1  real-time  polymerase  chain  reaction. 
DAC  =  2-deoxy-5-azacytidine. 


in  T  ceils  that  involve  alteration  in  methylation  of  isolated 
CpG  dinucleotides  in  gene-promoter  regions  devoid  of  CpG 
islands  [15,16],  Here  we  show  relative  demethylation  of  the 
PRV-1  gene  in  neutrophilic  granulocytes.  The  PRV-1  gene 
lacks  a  bona-fide  CpG  island;  there  are  only  14  CpG  sites 
in  the  1000-bp  region  upstream  and  200-bp  downstream 
of  the  transcription  start  site  (TSS,  Fig.  2).  Major  TSS  for 
the  gene  is  flanked  by  three  CpG  sites  located  30,  41,  and 
43  downstream  of  TSS.  We,  therefore,  analyzed  these 
CpG  sites  for  potential  relationship  with  gene  expression. 
Three  reported  alternate  potential  TSS  (UCSC  Genome 
Browser  on  Human  March  2006  Assembly)  are  positioned 
6-bp  upstream  or  8-bp  and  35-bp  downstream  of  the 
AF146747  transcript,  however,  the  three  studied  CpG  sites 
are  closest  to  those  TSS  as  well.  We  studied  methylation 
status  of  cytosines  30,  41,  and  43.  We  noticed  that  the  meth¬ 
ylation  status  of  C43  was  closely  correlated  with  the  PRV-1 
C42G  polymorphism  [13].  C43  methylation  levels  were  the 
highest  in  42G/G,  lower  in  42G/C,  and  the  lowest  in  42C/C 
genotype  (p  <  0.000 L,  Chi-square  test).  It  is  interesting 
that  42C/C  genotype  was  also  found  to  be  associated  with 
a  higher  expression  level  of  PRV-1  on  neutrophils  of  normal 
donors  [13].  One  might  explain  the  higher  expression  level 
of  PRV-1  in  42C/C  homozygotes  by  the  absence  of  the  C41 
CpG  site  and  by  a  decrease  in  methylation  of  C43  close  to 
the  initiation  codon  of  the  PRV-1  gene.  Whether  this  single 
nucleotide  polymorphism  affects  gene  expression  through 
alteration  in  the  DNA  methylation  status  is  unknown  and 
needs  additional  research,  however,  in  order  to  avoid  the  ef¬ 
fect  of  the  G42C  polymorphism  on  our  study,  we  focused 
our  additional  analysis  on  the  methylation  status  of  C30. 
We  found  that  C30  was  methylated  to  variable  extent  in 
all  DNA  samples  obtained  from  neutrophils.  Comparison 
of  PRV- 1-positive  neutrophils  to  PRV-1 -negative  neutro¬ 
phils  in  the  same  individual  showed  that  expression  of 
PRV-1  was  associated  with  a  lower  C30  methylation  of 
this  gene.  Cytosine  methylation  is  a  covalent  modification 
that  used  to  be  considered  as  a  DNA  alteration  with 
a  long-lasting  (and  perhaps  permanent)  regulatory  effect 
on  gene  transcription.  However,  there  are  examples  of  rapid 
demethylation  of  CpG  dinucleotides  in  a  gene  devoid  of 
CpG  island  at  TSS  after  activation  of  the  cell  by  proper 
stimuli,  such  as  rapid  demethylation  of  CpG  dinucleotides 
in  the  promoter  of  the  interleukin-2  gene  shortly  after  stim¬ 
ulation  of  T  lymphocytes  [17].  This  demethylation  is  inde¬ 
pendent  of  DNA  synthesis  and  increases  expression  of  the 
interleukin-2  gene  in  the  activated  lymphocytes.  Demethy¬ 
lation  of  the  PRV-1  gene  might  control  expression  of  this 
gene  under  physiologic  condition,  however,  exploration  of 
a  relationship  between  an  alteration  in  the  methylation 
status  and  overexpression  of  the  PRV-1  gene  detected  in 
patients  with  myeloproliferative  disorders  is  of  particular 
interest.  Overexpression  of  PRV-1  is  recognized  as  a  com¬ 
mon  event  in  patients  with  Philadelphia  chromosome-neg¬ 
ative  myeloproliferative  disorders,  and  especially  in  PV  and 
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ET.  We  found  a  relationship  between  the  mRNA  level  of 
PRV-1  and  the  methylation  status  of  its  gene  in  PV  and 
ET.  We  propose  that  methylation  of  the  PRV-1  gene  regu¬ 
lates  expression  of  this  gene  under  physiologic  conditions 
and  a  dysregulation  of  the  methylation  status  of  this  gene 
has  a  role  in  overexpression  of  PRV-1  in  PV.  Alternatively, 
demethylation  can  be  a  “passive”  event  that  follows  gene 
activation  (and  perhaps  augments  it)  rather  than  being  the 
initial  cause.  In  order  to  investigate  the  relationship  be¬ 
tween  the  level  of  DNA  methylation  and  mRNA  expres¬ 
sion,  we  manipulated  DNA  methylation  of  the  PRV-1 
gene  in  human  leukemic  cell  lines,  using  a  demethylating 
agent  decitabine  and  followed  the  change  in  the  mRNA 
level  of  PRV-1,  We  found  that  decitabine  demethylated 
the  PRV-1  gene,  which  was  associated  with  a  sharp  increase 
in  the  concentration  of  PRV-1  mRNA.  Several  recent  stud¬ 
ies  showed  an  important  role  for  JAK2V6,7F  somatic  muta¬ 
tion  in  pathogenesis  of  PV  [18-22].  Is  there  a  connection 
between  presence  of  JAK2V617F  and  level  of  methylation 
of  the  PRV-1  gene?  PRV-1  mRNA  expression  positively 
correlates  with  the  dosage  of  the  mutant  JAK2  allele  [23]. 
We  previously  showed  that  introduction  of  V617F  mutant 
JAK2  leads  to  increased  expression  of  both  PRV-1  protein 
and  mRNA  in  murine  myeloid  cells  [24].  Here  we  report 
that  PRV-1  gene  methylation  was  lower  in  the  presence 
of  activating  JAK2V6I7F  mutation.  One  explanation  for 
these  findings  is  that  JAK2  signaling  affects  DNA  methyl¬ 
ation.  We  hypothesize  that  augmented  IAK2  signaling  as 
a  result  of  JAK2V6I7F  mutation  alters  physiologic  activation 
of  neutrophils;  demethylates  the  PRV-1  gene  and  increases 
PRV-1  expression.  Supporting  a  possible  role  for  JAIC2 
signaling  in  epigenetic  regulation  of  the  PRV-1  gene 
expression  are  the  results  of  two  recent  studies  that  revealed 
JAK/signal  transducer  and  activation  of  transcription 
pathway  modulates  gene  expression  through  epigenetic 
mechanisms.  Shi  et  al.  [25]  showed  that  in  a  Drosophila 
melanogaster,  hematopoietic  tumor  model,  overactivation 
of  JAK2  globally  disrupts  heterochromatic  gene  silencing, 
and  Zhang  et  al.  [26]  demonstrated  a  role  for  phosphory- 
lated  signal  transducer  and  activation  of  transcription  3  in 
epigenetic  modulation  of  T-lymphocytes  genes  through 
inducing  activity  of  DNA  methyltransferase  1  gene. 
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Supplementary  Table. 


PCR  on  bisulfite-treated  DNA 
Forward  primer 
Reverse  primer 
Pyrosequencing  primer 
PCR  on  genomic  DNA 
Forward  primer 
Reverse  primer 
Pyrosequencing  primer 


5'-TGTTTAAGGGTTGGTATAAAG-3' 

5'-biotin-AACAAAATTCCCTTACCC-3' 

5'-AAAGAGATTATTAGTTATAG-3' 

5'-AGATTACCAGCCACAGACG-3' 

5'-biotin-GGAGGGCCAGCAGTAATA-3' 

5'-CAGACGGGTCATGAG-3' 


PCR  = 


polymerase  chain  reaction. 
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Abstract 

Decitabine’s  mechanism  of  action  in  chronic  myelomonocytic  leukemia  (CMML)  remains 
incompletely  understood.  We  studied  the  dynamics  of  neoplastic  cell  clearance  during 
decitabine  treatment  (1 00  mg/m2/course  every  4  weeks)  using  quantitative  monitoring  of 
mutant  alleles  by  pyrosequencing.  CMML  patients  were  first  screened  for  JAK2  and  NPM1 
mutations,  and  three  patients  with  mutations  were  identified.  Mutant  allele  percentages  in 
mononuclear  cell  DNA  were  followed  after  treatment,  along  with  methylation  of  LINE1  and 
ten  other  genes.  The  clearance  of  mutant  alleles  was  modest  after  the  first  cycle,  despite 
induction  of  hypomethylation.  Delayed  substantial  clearance  was  observed  after  2-4  cycles 
that  correlated  with  clinical  response.  Two  patients  had  complete  disappearance  of  mutant 
alleles  and  sustained  clinical  remissions.  In  another  patient,  mutant  allele  was  detectable 
at  clinical  remission,  which  lasted  for  8  months.  Our  data  suggest  a  predominantly  non- 
cytotoxic  mechanism  of  action  for  decitabine,  leading  to  altered  biology  of  the  neoplastic 
clone  and/or  normal  cells.  (1 50  words)  This  trial  was  registered  at  www.ClinicalT rials. gov 


as  #NCT00067808. 
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introduction 

Decitabine  has  therapeutic  activity  in  patients  with  myelodyspiastic  syndrome  (MDS)  and 
chronic  myelomonocytic  leukemia  (CMML).1,2  While  at  least  part  of  the  mechanism  of 
action  is  hypomethylation,  cytotoxicity  is  also  believed  to  play  a  role.  The  dynamics  of 
neopiastic  cells  during  treatment  may  be  assessed  with  appropriate  monitoring  of  genetic 
abnormalities.  Wijermans  et  al.  have  analyzed  the  dynamics  of  MDS  celis  after  decitabine 
treatment  by  following  cytogenetic  abnormalities.3  However,  this  method  requires  a 
baseline  cytogenetic  abnormality  and  multiple  bone  marrow  aspirations.  Furthermore, 
cytogenetic  analyses  typically  examine  only  20  celis,  which  is  not  suitable  for  sensitive 
quantification.  Molecular  genetic  abnormalities  can  become  useful  monitoring  tools  to 
overcome  this  issue.  Pyrosequencing  is  a  simple  method  to  quantitatively  detect 
nucleotide  polymorphisms.  The  benefits  of  this  method  are  simplicity  and  reproducibility. 
Here  we  studied  the  molecular  dynamics  of  CMML  in  three  patients  with  mutations  in  JAK2 
or  NPM1  who  were  treated  with  decitabine.  The  percentage  of  the  mutant  aliele  was 
quantified  and  stood  as  a  molecular  marker  of  disease  response  to  decitabine. 

Materials  and  Methods 

Peripheral  blood  samples  were  collected  after  informed  consent  from  16  patients  with 
CMML  upon  entry  to  a  phase  II  decitabine  study,  where  patients  were  randomized  to  one 
of  three  decitabine  schedules:  (1 )  20  mg/m2  intravenously  daily  for  five  days;  (2)  20  mg/m2 
subcutaneously  daily  for  five  days;  and  (3)  10  mg/m2  intravenously  daily  for  ten  days.2  The 
treatment  was  planned  to  be  repeated  every  four  weeks.  Blood  collection  was  scheduled 
on  day  0,  5,  12  and  28  during  the  first  cycle,  and  on  day  0  for  the  following  cycles.  DNA 
was  isolated  from  blood  samples  after  density  gradient  separation  of  mononuclear  cells 
using  standard  phenol-chloroform  extraction  methods.  Pretreatment  samples  were 
screened  for  mutations  of  JAK2  and  NPM1  using  pyrosequencing  as  previously 
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reported.4,5  This  research  was  approved  by  the  M.D.  Anderson  Cancer  Center  institutional 
Review  Board. 

JAK2  V617F  mutation  was  found  in  two  patients,  and  NPM1  mutation  was  detected 
in  another  patient.  In  patients  with  mutations,  the  proportion  of  mutant  alleles  was 
quantitatively  determined  by  pyrosequencing  using  samples  obtained  during  and  after 
decitabine  therapy.  Promoter  methylation  status  often  specific  genes  ( C1ORF102 ,  CDH1 , 
CDH13,  CDKN2B,  ESR1,  NPM2 ,  OLIG2,  PDLIM4 ,  PGRA,  and  PGRBf  that  are  frequently 
methylated  in  myeloid  malignancies  was  also  screened  in  pretreatment  samples  of  these 
three  patients  using  pyrosequencing.7  Methylation  status  of  the  LINE1  repetitive  element 2 
and  of  the  genes  that  showed  increased  methylation  before  treatment  was  followed  during 
and  after  treatment. 

Results  and  Discussion 

Patient  1  was  a  63-year-old  female  with  CMML  with  normal  karyotype  and  JAK2  mutation 
(1849G>T).  At  diagnosis,  white  blood  count  (WBC)  was  66  x  109/L  with  6%  monocytes  and 
4%  peripheral  blast  cells.  JAK2  mutation  was  detected  in  49%  of  the  alleles,  suggesting  a 
heterozygous  mutation.  After  the  first  cycle  of  decitabine  (schedule  2),  JAK2  mutant  alleles 
decreased  modestly  to  40%;  LINE1  analysis  showed  demethylation  from  59%  at  baseline 
to  49%  on  day  12,  and  then  re-methylation  to  56%  on  day  30.  After  the  second  cycle,  the 
patient  achieved  a  complete  hematologic  remission,  while  18%  of  alleles  still  carried  the 
JAK2  mutation.  Three  additional  courses  of  decitabine  cleared  the  mutant  down  to 
undetectable  levels  (<5%).  With  regards  to  methylation  changes,  CDH13  was  found 
methylated  (33%)  before  treatment.  At  clinical  complete  response  after  two  cycles  of 
decitabine,  CDH13  methylation  was  still  21%.  Three  additional  cycles  of  treatment 
decreased  methylation  down  to  the  level  seen  in  normal  controls  (10%)  (Figure  1).  The 
platelet  count  deteriorated  after  six  cycles  of  decitabine  despite  the  clearance  of  JAK2 
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mutant  alleles,  and  the  patient  was  taken  off  the  study  after  receiving  eight  cycles  of 
treatment.  The  patient  eventually  died  from  pneumonia  1 .4  years  after  initiation  of  therapy. 

Patient  2  was  a  78-year-old  male  with  CMML  with  normal  karyotype  and  NPM1 
mutation  (960-961  insTCTG).  At  diagnosis,  WBC  count  was  15  x  109/l,  with  13% 
monocytes.  He  had  anemia  and  thrombocytopenia.  NPM1  mutation  was  detected  in  35% 
of  the  alleles.  The  first  two  cycles  of  decitabine  (schedule  3)  were  associated  with  the 
typical  hypomethylation  induction  but  with  a  minimal  change  in  the  mutant  allele 
percentage  (from  35%  before  treatment  to  38%  after  two  cycles).  This  patient  achieved  a 
complete  hematologic  response  after  four  courses  of  decitabine  coincident  with  marked 
clearance  of  mutant  alleles  to  barely  detectable  level  (5%).  After  two  additional  courses  of 
decitabine  the  mutant  allele  was  no  longer  detectable.  PGRB  was  found  methylated  at 
baseline  (25%).  Although  transient  demethylation  of  PGRB  occurred  after  the  first  cycle, 
remethylation  was  observed  at  day  0  of  the  next  cycle.  Overall,  methylation  of  PGRB  in 
this  78-year-old  patient  followed  the  dynamics  of  LINE1  methylation  suggesting  that  PGRB 
was  moderately  methylated  also  in  normal  cells.  The  patient  remains  in  complete 
remission  at  2.8  years  after  initiation  of  therapy. 

Patient  3  was  a  55-year-old  male  with  CMML  with  normal  karyotype  and  JAK2 
mutation  (1849G>T).  At  diagnosis,  WBC  was  7x109/L,  with  16%  monocytes.  JAK2 
mutation  was  detected  in  35%  of  the  alleles.  LINE1  hypomethylation  was  induced  after  the 
first  cycle  (schedule  2)  with  transient  decrease  of  JAK2  mutant  alleles  to  20%.  The  second 
cycle  of  decitabine  was  associated  with  marked  clearance  of  the  mutant  cells  (6%  mutant 
alleles)  and  the  patient  achieved  CR  evidenced  by  disappearance  of  monocytosis  and 
normal  platelet  count.  Subsequently,  the  mutant  JAK2  allele  started  rising  again  (Figure  1), 
though  the  patient  clinically  remained  in  CR.  The  disease  overtly  progressed  after  8 
courses;  the  spleen  acutely  enlarged  and  bone  marrow  blast  count  increased  to  6%.  Given 
otherwise  stable  condition,  this  patient  received  one  more  course  of  decitabine  and  four 
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weeks  later  underwent  splenectomy,  which  confirmed  the  presence  of  CMML  in  the  spleen. 
The  patient  left  our  institution  after  splenectomy  to  receive  supportive  care  at  local  hospital. 
CDH1,  CDH13,  and  NPM2  were  found  methylated  at  baseline.  CDH1  and  NPM2  showed 
demethylation  after  three  courses  (from  17%  to  5%  and  from  28%  to  6%,  respectively), 
and  the  level  of  methylation  remained  low  until  the  obvious  progression  of  the  disease, 
except  for  one  point  of  transiently  increased  NPM2  methylation.  The  degree  of  CDH13 
methylation  was  20%  at  baseline.  Its  methylation  dynamics  essentially  followed  changes  of 
JAK2  mutant  allele,  suggesting  the  presence  of  methylation  in  malignant  cells. 

in  summary,  we  showed  that  clearance  of  neoplastic  cells  after  decitabine  therapy  in 
CMML  was  very  modest  after  one  cycle  despite  induction  of  LINE1  hypomethylation. 
Subsequent  cycles  were  associated  with  similar  LINE1  hypomethylation  dynamics,  but  a 
marked  delayed  clearance  of  the  mutant  clones  was  observed,  coincident  with  clinical 
remission.  Thus,  global  hypomethylation  precedes  clonal  elimination  and  clinical 
responses.  The  paucity  of  genes  hypermethylated  in  these  cases  makes  it  difficult  to 
comment  on  the  importance  of  tumor-suppressor  gene  demethylation.  Nevertheless, 
demethylation  dynamics  of  CDH13  gene  temporally  coincided  with  mutant  clone 
elimination.  Whole  genome  analysis  of  DNA  methylation  by  microarrays  or  deep 
sequencing  is  needed  to  uncover  more  gene-specific  methylation  changes  occurring  after 
decitabine  treatment.  Our  data  are  most  consistent  with  an  initial  modest  cytotoxic  effect  of 
decitabine  followed  by  a  delayed  and  complete  clearance  of  the  clone,  the  dynamics  of 
which  suggests  a  non-cytotoxic  mechanism.  Possibilities  for  this  delayed  action  include 
altered  biology  of  the  neoplastic  clone  (senescence  induction,  effects  on  stem  cell 
renewal),  induction  of  an  immune  response  8  against  CMML  clone,  or  effects  on  normal 
stem  cells.  Further  investigations  should  focus  on  these  mechanisms  in  order  to  improve 
the  results  seen  with  decitabine. 
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Figure  1  Legend 

Dynamics  of  DNA  methyiation,  mutant  alleles  as  neoplastic  clone  markers,  and 
peripheral  blood  counts. 

We  analyzed  methyiation  of  10  genes  at  baseline  in  all  three  patients,  and  genes  with 
increased  baseline  methyiation  were  followed  during  and  after  the  treatment.  Patient  1: 
CMML  with  JAK2  mutation  of  49%  alleles  prior  to  treatment,  which  decreased  to 
undetectable  level  after  4  cycles  of  decitabine.  CDH13  methyiation  was  present  at  initiation 
of  therapy,  and  decreased  after  decitabine.  Patient  2:  CMML  with  NPM1  mutation  of  35% 
alleles  prior  to  treatment,  which  decreased  to  undetectable  levels  after  4  cycles.  PGRB 
methyiation  followed  the  dynamics  of  LINE1  methyiation  after  decitabine.  Patient  3:  CMML 
with  JAK2  mutation  of  35%  alleles  prior  to  treatment,  which  decreased  to  6%  after  2  cycles 
of  decitabine  and  then  increased  again.  CDH1,  CDH13  and  NPM2  methyiation  was 
present  at  initiation  of  therapy,  and  decreased  after  decitabine.  Cycles  of  decitabine 
treatment  are  indicated  as  purple  rectangles  above  the  graphs. 


